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CHAPTER  1 


INTRODUCTION 

In  numerous  applications,  there  is  a pressing  need  for  a 
short  distance  HF  antenna  which  exhibits  desired  electrical 
characteristics  but  which  also  possesses  acceptable  mechanical 
features.  For  an  HF  antenna  to  be  of  greatest  utility,  its  input 
impedance  and  VSWR  must  not  vary  outside  of  prescribed  bounds 
over  the  frequency  band,  and  its  radiation  pattern  must  be 
essentially  omni-directional  throughout  this  frequency  range 
of  interest.  The  antenna  must  be  quickly  erectable,  highly 
portable,  light  weight,  and  must  be  able  to  withstand  natural 
environmental  hazards  such  as  wind,  ice  loading,  rain,  corrosion, 
etc.  Finally,  the  antenna  must  be  simple,  of  low  cost,  and  its 
properties  under  varying  ambient  conditions  must  be  well  known. 

Of  course,  to  achieve  such  an  ideal  radiator  is  impossible, 
but  it  certainly  is  desirable  to  attempt  to  realize  a physical 
structure  having  as  many  of  the  above  features  as  is  practical. 

It  is  apparent  that  a careful  analysis  must  precede  the 
development  of  an  optimum  HF  antenna,  since  any  design  task 
necessitates  intelligent  compromises  involving  trade-offs  and 
modifications  among  ideal  characteristics.  Only  after  full 
understanding  of  a generic  antenna  structure  and  a delineation 
of  its  properties  and  their  interrelationships  can  an  engineer 
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undertake  the  difficult  task  of  designing  a structure  which 
is  near  optimum  for  his  applications. 

A basic  radiating  element  which  can  be  made  to  operate 
satisfactorily  over  a broad  range  of  frequencies  is  the  so- 
called  biconical  antenna. 

Biconical  dipoles  have  been  used  extensively  as  broad- 
band antennas  in  the  VHF-UHF  frequency  range.  The  biconical 
type  structure  in  free  space  has  been  studied  by  various 
authors  as  a spherical  b ound a r y- val ue  problem,  and,  from  this 
point  of  view,  input  impedance  can  be  predicted  from  analytic 
formulas  for  small  cone  angles.  Impedance  data  for  large 
angle  biconical  antennas  is  derived  mainly  from  experimental 
studies,  and  it  is  well-known  from  these  results  that  the  input 
resistance  and  reactance  remain  fairly  constant  over  a wide 
frequency  range  for  cone  angles  from  20°  to  90°as  compared  to 
the  thin  cylindrical  dipole.  The  gain  of  the  biconical  struc- 
ture also  exhibits  this  resonably  constant  behavior  as  a 
function  of  frequency.  Experimentalists  have  found  that,  by 
adjusting  the  total  length  and  the  angle  of  the  conical 
sections,  optimum  radiation  characteristics  can  be  obtained 
for  frequency  range  variations  of  several  orders  of  magnitude. 

In  the  size  needed  for  operation  at  HF,  the  solid  or 
shell  biconical  antenna  would  be  so  massive  that  it  would  be 
totally  impractical.  However,  the  properties  of  the  bicone 
are  so  near  ideal  compared  with  those  of  othersimple  antennas, 
one  is  naturally  led  to  investigate  the  possibility  of  developing 


a biconical-like  antenna  which  retains  desired  electrical 
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features  but  which  is  mechanically  realistic.  A simple  radi- 
ator which  may  be  viewed  as  a form  of  a blconical  antenna  can 
be  fabricated  from  conducting  wires  having  the  configuration 
shown  in  Fig.  1-1.  More  sophisticated  antennas,  comprising 
various  combinations  or  arrays  of  elemental  wire  bicone  radi- 
ators, may  be  desired  where  prescribed  applications  warrant 
added  complexity. 

Even  though  the  classical  biconical  antenna  is  known  to 
exhibit  a remarkably  wide  bandwidth,  how  similar  are  the 
properties  of  the  wire  structure  depicted  in  Fig.  1-1  to 
those  of  the  classical  radiator?  Does  it  exhibit  broadband 

This  report  describes  a research  program  in  which  are 
developed  the  tools  for  investigating  the  wire  biconical 
antenna  of  Fig.  1-1  in  free  space  as  well  as  above  an  ideal 
ground  (three  orientations).  A new  theory  for  analyzing  wire 
structures  has  been  developed  and  is  described  in  this  report. 
A numerical  solution  procedure,  based  upon  the  new  theory 
and  applicable  to  general  wire  structures,  is  discussed  and 
sample  calculated  data  are  presented.  A computer  program 
has  been  developed  for  implementing  the  theory  and  numerical 
solution  procedure.  Also,  calculations  have  been  made  for 
various  configurations  of  the  wire  biconical  antennas,  and 
selected  results  are  presented  and  compared  with  available 
measured  data. 
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CHAPTER  2 


FORMULATION  OF  STRAIGHT-WI RE  INTEGRAL  EQUATIONS:  PIECEWISE 

LINEAR  AND  PIECEWISE  SINUSOIDAL  TESTING 


2.1  INTRODUCTION 


In  this  section  a particular  moment  method  [1]  analysis 
of  s t rai ght -wi re  scatterers  and  antennas  is  developed,  and  the 
associated  integral  equation  for  the  straight  element  is  form- 
ulated. The  method  emphasizes  piecewise  linear  and  sinusoidal 
testing,  and  interrelationships  between  these  two  testing  pro- 
cedures are  delineated.  Both  pulses  and  triangles  [1]  are  used 
for  representation  of  the  current  and  a brief  discussion  of  the 
utility  and  advantages  of  each  is  provided.  Example  results  of 
straight-wire  calculations  are  included. 

From  basic  electromagnetic  theory  applied  to  the  perfectly 
conducting  cylinder  of  Fig.  2.1, assumed  to  be  excited  by  some 
incident  field,  one  may  readily  derive  the  following  fundamental 
i n t e g r o- d i f f e r en t i al  equation. 


( -^  + r2)  a^(z)  = -2^  E^(z)  , 


(2.1) 


relating  the  magnetic  vector  potential, 

L/2 


= A / i 


i(z')K(z-z')dz'  , 


z L/2 


(2.2) 


’•srx-rr 


Fig.  2.1.  Straight-Wire  Scatterer  Illuminated  by 
Incidemt  Field 


Fig.  2.2.  Straight-Wire  Antenna  Excited  by  Voltage 
Source  Located  at  z = z g 


h 


.jrf.  'r  -■  ' 


Co  the  unknown  incident  electric  field  having  an  axial  com- 


ponent on  the  surface  of  the  cylinder,  which  resides  in  a 
homogeneous  space  characterized  by  (y,e,o=0).  The  kernel  in 
(2.2)  ab  ove  i s 


K(0 


+ sin^4’’/2]^ 

e ^ ^ 

+4a  sin  4>'/2]  ^ 


(2.3) 


where  k is  27t/A  (wavelength)  at  the  angular  frequency  to  of  the 
suppressed  harmonic  time  variation  Through  (2.2),  Equa- 

tion (2.1)  relates  the  known  excitation  to  the  unknown  axial 

z 

current  i on  the  scatterer  which  is  specified  to  be  a tube  of 
radius  a and  length  L as  suggested  in  Fig.  2.1.  If  one  wishes 
to  treat  the  antenna  problem  in  which  the  excitation  is  a slice 
generator  of  voltage  V impressed  at  a point  z , he  replaces  E^ 
of  (2.1)  by  V6(z-z  ) where  5(z)  is  the  familiar  delta  function 
(Fig.  2.2). 

For  present  purposes,  the  wire  radius  is  looked  upon  as 
being  very  small  relative  to  the  wavelength  A as  well  as  to  the 
cylinder  length.  Such  restrictions,  common  in  thin-wire  analyses, 
assure  one  that  the  current  on  the  cylinder  is  circumferentially 
independent  and  that  it  can  be  accounted  for  by  the  total  axial 
current  i.  Under  these  thin-wire  assumptions,  one  could  make 
use  of  the  so-called  reduced  kernel  approximation  to  (2.3), 
but,  since  the  solution  procedures  below  are  not  unduly  com- 
plicated by  the  exact  kernel,  K(0  of  (2.3)  is  used  without 
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approximation.  Even  though  thin  wires  are  under  investigation 


here  the  analysis  applies  to  a thick  cylinder,  if  only  an  axially 


directed  current  exists  and  is  circumferentially  independent. 


as  would  be  the  case  for  an  Isolated  antenna  driven  by  a cir- 


cularly uniform  source  or  generator. 


The  solution  procedure  presented  here  is  essentially  the 


moment  method  [ 1],  but  a change  in  the  sequence  of  the  method’s 


usual  steps  enables  one  to  gain  interesting  insight  into  the 


nature  of  the  numerical  procedure.  The  alteration  alluded  to 


is  that  one  performs  the  testing  of  the  equation  to  be  solved 


(2.1)with  elements  T^  of  the  selected  testing  set  before 


the  unknown  current  is  approximated  as  a linear  combination  of 


the  elements  1 of  the  basis  set  ( i 'S 
n \ n J 


2.2  TESTING 


To  ensure  the  stability  of  the  anticipated  system  of 


linear  equations  to  be  formed  in  the  solution  procedure,  only 


subdomain  testing  sets  are  considered  here.  Pursuant  to  such 


a system,  one  equates  the  corresponding  projections  of  the  two 


sides  of  (2.1)  onto  the  space  spanned  by  , which  implies 


that  the  scalar  products  formed  by  T^  and  the  two  sides  of 


(2.1)  are  equated  for  each  m; 


J — E , T 
ui  z m 


(2.4a) 


.J 


wliere  tlie  scalar  product  is 


. g>  = 


f g*dz 


(2.  Ab) 


For  real  testing  sets  subdomain  type  considered 

here,  where  T (z)  is  zero  outside  ze(z  .),  (2.4)  is  simply 

X!i  TI\—  i iut  i 

^m+l  ^m+1 


(2.5) 


The  two  testing  sets  which  are  particularly  amenable  to  num- 
erical solutions  of  (2.1)  are  the  piecewise  linear  and  piecewise 


^ o ■Tocnrfanf'ivplv. 


where  2A  = (z  , - z , The  former  set  above  is  illustrated  in 

m+1  m-1 

Fig.  2.3  where  one  observes  overlapping  triangles  on  the  interior 
intervals  with  h a 1 f - 1 r 1 ang le s (dashed)  on  the  lower  and  upper 
end  intervals.  Testing  with  the  end  h a 1 f - 1 r 1 an g le s , or  end 
"half-sinusoida"  in  the  case  of  piecewise  sinusoidal  testing,  in 
the  moment  method  admits  the  representation  of  the  current  in 
a series  which  does  not  satisfy  a priori  the  boundary  conditions, 


i(+L/2)  = 0 


(2.7a) 


i(-L/2)  = 0 


(2.7b) 


and  provides  a convenient  way  to  incorporate  the  boundary 
Conditions  (2.7)  of  the  problem.  Having  a means  to  Incorporate 
conditions  at  wire  ends  becomes  important  in  situations  where 
wires  join. 

Interior  Testing  Testing  of  (2,1)  with  the  elements  (2.6a) 

of  the  piecewise  linear  set  over  the  interior  full  intervals 

of  length  2A,  (z  ,,z  .,),  m = 0 ± 1,  ± 2 ,..,+[ M- 1 ] , yields 
m-i  m+i 


(See  Appendix  I ) 


'm+1 


A (z  ,)  - 2A  (z  ) + A (z  ,)1  + 1 A_(z)A^(z)dz 

z tn  I X z tn  z ni”  X j ■ 


z ra  z 


z m 


Z*2 


m-1 


(2.8a) 


.k‘ 


= -J- 


m+1 

/ 


(z)A  (z)dz 
m 


z = z 


m-1 
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and  (2.'J\>)  becoue,  respectively. 


'-M+  1 


d z z - M 


. k‘ 


= -J- 


and 


/ 


/.\(z)/. 

-M  Z 


■-M 


-M+  1 

/ 


, 1 


■ „ ( 2 ) E 
-M  z 


d z z M 


r 


'M-1 


(z ) (z ) dz 


Also,  for  end  testing  with  hal f-s inusoids , one  obtains 


d , / , , cos  IcA  . , , k 

dz‘  z ^-M  sin  kA  ^ s in  kA  ^z^^-M+1^ 


= 'J 


0) 


/ 


■-M+1 


z = z 


(z  ) A_j^  (z  ) dz 


-M 
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( Z ) dr 


( z ) dz 

(2  . 10a) 


. (2 . 10b) 


(2.11a) 


4 

\ 


i| 

i 


and 


sin 


kA 

kA 


A..  (Zm)  + 


77a  a (z„_.) 


( z ) ( z ) dz  , (2.11b) 

2 M 

M-  1 

^ « s 

where  T ^ may  be  either  or  . 

2.3  BOUNDARY  CONDITIONS 

Attention  is  called  to  the  derivatives  of  the  vector 
potential  at  z + which  appear  in  Eqs.  (2.9)  and  (2.10). 

Appealing  to  the  Lorentz  condition,  one  can  readily  relate 
these  terms  to  the  electric  scalar  potential  ♦ by  use  of 


_d_ 
d z 


A 

z 


(2.12) 


I 


I 


) 

where,  of  course,  <t  is  calculated  from  the  electric  charge  on 
the  cylinder.  If  one  does  not  enforce  i(z^j^)=0  a priori,  then, 
• in  order  that  the  charge  and  current  on  the  cylinder  constitute 

1 a complete  electric  system  [2],  rings  of  discrete  charge  must 

reside  on  the  upper  and  lower  rims  of  the  tube.  Since  the 
^ current  is  assumed  uniform  over  any  cylinder  circumference, 

i 

1? 

« 
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f . i c li  r i n f li  .1  r is  u n 1 1 o r m an  cl  has  density* 


j 2VUiB 


where  the  upper  and  lower  signs  are  associated  with  quantities 
at  and  z_^  (L/2  and  -L/2),  respectively.  Now  one  may  write 

the  expression  below  for  the  scalar  potential  at  a point 
z£;(-L/2,  l/2)  on  the  cylinder  in  the  form 


^ f 


— i (z')  K(z-z')dz' 


where  the  integral  contribution  to  the  potential  Is  due  to 
the  distributed  charge  and  the  remaining  two  terms  come  from 


the  ring  charges.  The  potential  is  to  be  evaluated  on  the 
cylinder  at  at  which  points  (2.13)  is  unbounded.  Specifically 

one  can  show*that  the  potential  exhibits  the  following  properites 


at  the  upper  and  lower  ends,  respectively: 


+ 

1>  = liim  'i’Cz)  = 2 

z f Z-.  i 4tt  coe a z +z.. 

M M 


Ht!)]' 


(2 . lAa) 


*See  Appendix  II. 


m 


and 


= ilim  <{>(z) 


2 

jAiT  wea 


£ 1 m 
*-  z 4-2 


£n 


2a 


U-2-m 


ll 


-M 


(2.14b) 

Even  though  (2.14a)  and  (2.14b)  are  unbounded,  the  derivatives 
of  vector  potential  in  (2.10)  and  (2.11)  are  replaced  by 
equivalent  terms  involving  explicitly  in  order  to  facilitate 

discussions  which  follow. 

2.4  approximations 


For  a function  f sufficiently  smooth  over  the  interval 

( 2 z ) one  employs  the  approximation, 

m- 1 ’ ro+1 


'm+1 


/ 


f(z)T  (z)dz  = Af(z  ) 
m m 


(2.15) 


z-z 


in~  1 


i 

V 


9^  s 

where  T here  may  be  either  A or  A . In  a numerical  solution 
m mm 

procedure,  A must  be  small  compared  to  the  wavelength  A,  if 
(1)  the  excitation  is  to  be  sampled  adequately  along  the  wire 
and  (2)  the  number  of  elements  in  the  basis  set  {‘4  for  the 
current  approximation  (hence  the  number  of  equations  obtained 
from  testing  Equation  (2.1))  is  to  be  large  enough  to  accurately 


lb 


J 


A similar  approximation  may  be  made  for  Lestinji 


r !•  p r f s e n t i . 

over  the  end  intervals  with  the  right  sice  of  (2.15)  replaced 
by  since,  for  end  testing,  the  integration  is  over 

an  interval  of  length  A rather  than  2A  as  is  the  case  in  (2.15). 


2.5  SUMMARY  OF  TKSTING 

Employing  the  approximations  above  in  (2.8)-(2.11)  and 
making  use  of  (2.12)  and  (2.14)  in  (2.10)  and  (2.11),  one 
arrives  at  the  two  systems  of  equations  below  for  piecewise 
linear  and  piecewise  sinusoidal  testing.* 


Piecewise  Linear  Testing 


1- 


(kA)' 


A ^z  [M-1]  ^ 


.k^  A pi. 

2 ’^z^"- 


(2 . 16a) 


A 

m = 0,  ±1,  +2, i[M-l] , 

(2. 16b) 


(kA)' 


I A (z  ) + A (z  ,)  = -j — AE^(z  ) 

f z m z'm-lj  z m 


1- 


(kA)' 


*Another  form  of  the  equations  resulting  from  end  testing 
is  convenient  for  some  purposes;  this  form  is  discussed  in 
Chapter  3. 
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Pie  wise  ^ _i_n  u soidal  Test! 


, cos  kA  , , , . k , 

sin  kA  z^^-M  sin  kA  z ^-[M-1] 


V . . k A , , i , - 

^ = "^17  2 


(2.17a) 


sin  k A 


A (z  - 2cos  kA  A (z  ) + A (z 

z m+1  z m z ra 


-j~  AF>(z  ), 
u z m 


m = 0,  ±1,  ±2,. . . , ±[M-1] , 


(2.17b) 


, cos  kA  . ^ , 

k — : nr  A (z„ ) + 

s 1 n kA  z M 


sin  kA  ^z  [M-1 ] ^ 


• A „i  , 


(2.17c) 


2.6  RELATIONSHIPS  BETWEEN  PIECEWISE  LINEAR 
AND  PIECEWISE  SINUSOIDAL  TESTING 

For  sufficiently  small  A,  one  observes  that 

J?.  s 

and  A do  not  differ  significantly,  and,  therefore,  operations 
ra  ra 

with  one  testing  element  should  produce  results  similar  to 
those  achieved  via  the  other.  Such  is  indeed  the  case  as  is 
seen  from  the  discussion  below.  To  obtain  (2.16)  from  (2.8a) 
and  (2.10),  approximation  (2.15)  is  applied  to  two  Integrals 
in  each  equation,  whereas  (2.17)  follows  from  (2.8b)  and  (2.11) 


Miljjcit  to  approximation  of  a sin^lo  inte}>rai  per  equation. 


in  this  sense,  Eqs.  (2.17)  resulting  from  piecewise  sinusoidal 
testing  are  more  nearly  exact  tlian  are  (2.16)  resulting  from 
piecewise  linear  testing.  Of  course,  in  a numerical  solution 
one  may  make  use  of  Eqs.  (2.8),  (2.10),  and  (2.11)  which  are 

exact,  but,  as  one  would  suspect  in  view  of  the  closeness  of 

* S 

to  for  A<<>,  (2.16)  and  (2.17)  are  seen  to  be  the  same 

set  of  equations  within  the  approximations. 


i 

f 


sin  kA  i kA  , 


cos  kA 


(kA)^ 

2 


which  are  quite  good  whenever  kA<<l.  Of  additional  interest 
is  the  observation  that  (2.16b)  is  the  difference  equation 
approximation  to  (2.1)  arrived  at  by  replacing  the  second  deriv- 
ative operator  by  the  second  central  difference  operator.  Thus, 
piecewise  linear  testing  of  Eq . (2.1)  is  approximately  equivalent 

to  replacing  it  by  the  corresponding  difference  equation. 


2.7  USE  OF  PIECEWISE  LINEAR  AND  PIECEWISE 
SINUSOIDAL  TESTING  IN  NUMERICAL  SOLUTION  OF 
CURRENT  ON  STRAIGHT-WIRE  SCATTERER 

At  this  point  attention  is  given  to  the  task  of  solving 
for  the  unknown  current  i(z)  induced  on  a s t r ai gh t -w i r e scat- 
terer  by  a specified  incident  field  E^.  To  obtain  such  a 
solution  one  establishes  a procedure  based  upon  either  (2.16) 
or  (2.17).  In  the  procedure  the  unknown  current  i(z)  is 

20 


approximated  as  a linear  combination  of  known  elements  l^(z) 
of  a suitably  chosen  basis  set  (j- ^ • 


(2.18) 


where  the  I 's  are  unknown  complex  constants  to  be  calculated 
n 

in  the  procedure.  With  i(z)  of  (2.18),  the  vector  potential 
of  (2.2)  can  be  written 


A (z)  = 7 I A (z) 

z Z-  " ^n 


(2.19a) 


where 


A^  (z) 
n 


I i^(z')K(z-z')dz' 


(2 . 19b) 


is  a partial  vector  potential  contributed  by  the  basis  element 

i . Once  a basis  set  < i > has  been  selected  in  the  numerical 
n V n y 

solution  technique,  i^  would  be  known  and  one  could  immediately 

determine  A . Next  one  substitutes  (2.19)  into  the  tested 
z 

n 

i n t e g r o- d i f f e re n t i a 1 equation  ((2.16)  or  (2.17)).  For  illus- 
trative purposes  (2.16)  is  employed  here  and,  to  facilitate 
explanation,  one  finds  the  following  definitions  to  be  convenient 


i: 

V 

* 


L 


2(1-  iMi  )k(z^-z')  + K<z_ 


dz  ' 


(2 . 20a) 


and 


V = -j— AE^(z  ) . 

m u z in 


(2 . 20b) 


In  terms  of  the  above  definitions,  (2.16b)  may  be  written 
compactly  as 

N 


I 

n = -N 


IS  = V 
n mn  m 


m = 0,  ±1,  ±2,...,  ±(M-1)  , (2.21) 


which  is  a system  of  (2M-1)  linear  equations  with  (2N+1)  un- 
knowns (I  's).  To  satisfy  the  boundary  conditions, 
n 

i(-L/2)  = i(L/2)  = 0,  one  of  the  following  approaches  may  be 
emp loyed . 

1.  Equations  (2.16a)  and  (2.16c),  with  (2.18) 
and  (2.19),  may  be  appended  to  (2.21)  in 
which_^case  the  presence  of  the  terms  <5> 
and  $ of  (2.16)  in  such  an  augmented  system 
of  equations  forces  satisfaction  of  the 
boundary  conditions. 
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r 


2.  The  boundary  conditions  imply 

N 


i(-L/2)  = >I  i^(-L/2)  = 0 

n n 


n = -N 


and 


I 


i(L/2)  = 7 I^i^(L/2)  = 0 

n = -N 


which  two  equations  may  be  added  to  the 
sys  tern  (2.21) . 


3. 


If  the  basis  set 


V ny 


is  of  the  subdomain 


type  [1]  , the  boundary  conditions  may  be 
incorporated  directly  by  setting  i ^ and 

i-,  equal  to  zero;  i „ and  i., 

N -N  N 


exist  ove  r the 
lower-end  and  upper-end  subdomains,  respectfully, 


With  the  wire  boundary  conditions  honored,  the  system  of 
equations  (2.21)  can  be  solved  for  the  I^'s,  which  in  conjunc- 
tion with  (2.18)  lead  to  an  approximation  for  the  wire  current 
i (z)  . 

Pulse  Basis  Set  As  an  example  of  how  the  procedure  described 
here  may  be  implemented,  we  choose  pulses  [1]  for  the  elements 
of  the  basis  set,  in  terms  of  which  (2.18)  becomes 


i (z  ) = 


I 


n = -N 


^nPn^") 


(2.22) 


where  the  pulse  function  p^(z)  is  given  by 
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(2.23) 


r 1,  zt(i;^-A/2,  r^^^  + M2) 

P^(z)  = I 

/ 0,  z«'(c  -A/2,  C +A/2)  . 
^ n n 

and  where 


= nA  , n = 0,  ±1,  +2,...,  ±N  . (2.24) 


In  Fig.  2.4  a current  approximated  by  pulses  is  depicted  and 

C of  (2.23)  is  shown  as  the  coordinate  location  of  the  center 
n 

of  the  n*”^  pulse.  Notice  in  Fig.  2.4  that  and  1^^  are 

made  zero  in  order  that  boundary  conditions  are  satisfied  by 
the  current  approximation;  thus,  (2.22)  reduces  to 


N-1 

n=- (N-1) 


(2.25) 


and,  with  i replaced  by  p in  (2.20a),  S becomes 

Ti*  ''*r»  'mi-* 


S 

mn 


1 _H_ 
A 4tt 


z 


K ( z - z 
m 


) 


+ (2.26) 

since  p„(z)  is  zero  outside  ze(i;  -A/2,  Z,  +A/2).  A simple 
n n n 

change  of  variables,  use  of  (2.24)  , and  replacement  of  z by 

m 

(see  Fig.  2.4) 


2 4 


mA  , in 


0 , • J , '2  


(2.27) 


allow  one  to  convert  (2.26)  to 


A/2 


47t  A 


/ 


K (A [ra+l-n ] + C ' ) 


2(1-  1 K(A[m-n]  + C'  ) 


dC’ 


+ K ( A [ m- 1 -n ] ' ) 


(2 . 28) 


a form  which  Is  computationally  more  convenient  than  (2.26). 
With  N=M  and  with  (2.25)  and  (2.28)  in  (2.21),  one  has  a system 
of  (N-1)  equations  in  (N-1)  unknowns  from  which  he  can  deter- 
mine the  I 's  (unknowns).  ^'-g.  2.5  shows  a wire  of  length  L 

n 

partitioned  into  10  segments  of  length  A with  the  current 
approximated  by  9 pulses.  Notice  the  geometric  relationships 
between  the  location  of  the  piecewise  linear  testing  functions 


. i. 


and  of  the  pulses  and  observe  that  the  current  on  each  end 


is  represented  by  half-pulses,  each  of  which  is  made  zero  so 
tliat  the  current  approximation  satisfies  the  boundary  condition 
from  the  outset.  From  this  figure  we  generalize  as  follows. 

For  a wire  of  length  L one  may  select  any  number  P (odd  number 
for  convenience)  of  unknown  pulses  and  partition  L into  (P+1) 
segments  of  length  A=L/(P+1),  which  lead  one  to  specify  the 
number  of  equations  in  the  linear  system  (2.21)  to  be,  of  course, 
P , 

P - number  of  unknown  (pulses)  and  number 

of  equations  in  (2.21) 

A - L/(P+1) 
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Fig.  2.5.  Relative  Uxations  of  Pulses  (Bases  Set)  and  Triangles  (Testing  Set)  on  Straight  Wire 


.Numerical  Calculatlon.s  We  observe  from  (2.3)  and  (2.28)  that, 
when  m and  n are  related  by  either  m=n-l,  m=n,  or  m=n+l,  one 
of  the  tliree  integrands  of  (2.28)  is  singular  leaving  one  with 
the  need  to  evaluate  integrals  like 


K(C')dc’ 


i;'--A/2 


with  a singular  integrand.  This  integral  is  discussed  and 
evaluated  in  a forthcoming  paper  [3]  by  Butler.  All  com- 
binations of  n and  m other  than  those  cited  above  give  rise 
to  nonsingular  integrands  of  (2.28)  which  may  be  integrated 
numerically  with  great  ease  and  accuracy. 

Piecewise  Linear  Basis  Set  As  another  example, one  may  consider 
— - « 

representing  the  current  as  a linear  combination  of  piecewise 

a z 

linear  functions  , where  the  ' s of  (2.6a)  replace  the 

P^'s  of  (2.18).  These  are  often  called  triangles  or  triangular 
pulses  and  are  illustrated  in  Fig.  2.6  where  one  sees  the  cur- 
rent approximated  by  9 triangles  and  tlieir  relative  positions 
t o the  9 testing  triangles.  The  peak  of  a given  triangle 

used  in  the  s’t  to  approximate  current  is  equal  to  current  at 

S, 

the  point  where  the  peak  occurs,  i.e.,  the  peak  of  is  I^, 

and  triangles  are  seen  to  be  arranged  in  sucli  a way  that  the 
boundary  conditions  are  satisfied.  Also,  an  approximation  of 
current  by  a sum  of  weighted  triangles  is  equivalent  to  piece- 
wise  linear  interpolation  as  suggested  by  the  dashed  lines  of 
Fig . 2.6. 


£ 

The  coefficients  S are  needed  with  i =A  in  order  that 

mn  n n 

one  mav  establish  the  linear  system  (2.21)  when  the  current  is 

represented  by  a piecewise  linear  approximation.  To  obtain 

£ 

S with  i =A  , A of  (2.19b)  is  determined,  reduced  to  a use- 
mn  n n 2 

n 

ful  form,  and  then  substituted  into  (2.20a).  First,  from  (2.19b), 

I./2 


A (z  ) = 7^  (A  (z  )K(z  -z  )dz 
z m 4tt  I n m 

n J 


z '=-L/2 


and,  since  A is  nonzero  only  over  z'e(C  +A),  the  above 

n n n 


expression  becomes 


(n+l)A 


Z 4TT  ; n 

n 


/ 


A (mA)  = I A^  (z  ' )K(Am-z  ' )dz  ’ 


z ' = (n-1) A 

where  use  is  made  of  (2.24)  and  (2.27).  Subject  to  the  change 
of  variables  z ' = - c' +n A , on e has 


Li 


A (mA)  = 7^  /A^(nA-C')K(A[m-n]+^')dS’ 

z 4tt  j n 

n 

?’=-A 

and,  as  can  be  seen  from  (2.6a)  or  (1-1), 


A^;(nA-?')  = ^[A-U'  I ] 


= Aq(^') 


so  one  finally  arrives  at 
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A 


A ( mA  ) 
z 

n 


P 

4tt 


(C’)K(A[m-n]+^’)dc’ 


C'--A 


(2.29) 


Substitution  of  (2.29)  into  (2.20a)  enables  one 
for  a piecewise  linear  approximation  of  i(z): 


to  obtain  S 

mn 


A 

^mn  = ^ K(A[m+l-n]  + ^')  -2(1-  -^2—^)  K(A[m-n]+C’) 

?'=-A 

\ 

+ K(A[m-l-n]+C' ) J dC  . (2.30) 

Equation  (2.30)  is  in  a highly  convenient  form  for  numerical 
'i|  computation.  As  is  the  case  in  (2.28),  the  kernels  of  (2.30) 

f can  be  singular;  such  integrals  are  discussed  in  [3]  and  sugges- 

tions for  evaluation  are  given. 

^ With  available  for  the  current  approximated  by  triangles, 

. one  substitutes  (2.30)  into  (2.21)  and,  with  N=M,  solves  for  the 

^n'®'  Knowing  these  coefficients,  one  approximates  the  current 
by  means  of 


(N-1) 


t 


4 


i 

V 


Samplf  Kobalts  TIte  .ib(jve  theory  is  directly  applicable  to  the 

problem  ol  calculation  ol  current  on  a straight-wire  scatterer. 

By  solving  the  set  ol  algebraic  equations  (2.21)  with  given 

by  (2.2B),  one  obtains  the  l„'s  of  the  pulse  current  expansion 

(2.25)  and,  lienee,  has  an  approximation  for  the  current  i(z); 

similar! V,  solution  of  (2.21)  with  S of  (2.30)  yields  the 

mn 

I 's  associated  with  a triangle  expansion  which  may  be  used 
n 

in  (2.31)  for  determination  of  current.  As  an  example  of  re- 
sults obtained  from  the  present  analysis.  Fig.  2.7  depicts  the 
current  induced  by  a normally  incident  field  of  1 volt/meter  on 
a straight-wire  scatterer  of  length  l.=  A/  4 and  radius  a=A/1000. 
•Notice  that  data  obtained  by  pulse  and  triangle  current  expan- 
sion are  essentially  identical. 

Convergence  in  order  that  one  may  observe  the  convergence  of 
the  series  (2.25)  'and  (2.31)  for  s t r a i gh  t -w  i re  scatterers  of 
various  lengths,  all  of  radius  a=A/1000  and  illuminated  by  a 
normally  incident  electric  field  of  1 volt/meter,  Tables  I and 

II  are  provided.  From  a careful  study  of  these  tabulated  data, 
one  sees  little  difference  in  convergence  rates. 

Other  Basis  and  Testing  Sets  One  may  choose  piecewise  sinu- 
soidal testing  functions  and  employ  Equation  (2.17b)  rather 
than  (2.16b)  to  solve  for  the  current  on  the  scatterer.  A 

procedure  based  upon  (2.17b)  would  parallel  those  described 

s fL 

above.  However,  in  view  of  the  closeness  of  A to  A for  suf- 

m ra 

ficiently  small  A/A,  one  would  obtain  essentially  the  same 
numerical  results.  Furthermore,  one  may  wish  to  specify  the 
elements  of  the  basis  set  to  be  but,  here  again,  tlie  close- 

3 1 


Fig.  2.7.  CuRREiir  on  Straight-Wire  Scatterer 


BLE 


*Number  of  Subdomains  is  One  More  than  Number  of  Current  Elements 
**Absolute  Value  in  Parenthesis 


TABLE  II 


Y/1)  pazTXPmaoN 


of  Subdomains  is  One  More  than  Number  of  Current  Elements 


ness  of  lo  suggests  that  a piecewise  sinusoidal  approxi- 

mation of  current  would  lead  to  solutions  wliich  differ  only 
insignificantly  from  those  obtained  via  a piecewise  linear 
approximation.  In  summary,  then,  the  two  examples  given  above 
of  solution  procedures  can  be  viewed  as  covering  the  following 
combinations  of  t e s t i n g / c u r r e n t approximation:  (1)  piecewise 

linear /pulses , (2)  piecewise  1 in e a r / p i e cew i s e linear,  (3)  piece- 

wise  1 ine ar / p ie cewi s e sinusoidal,  (4)  piecewise  s inus oi dal / pulses  , 
(5)  piecewise  sinusoidai/piecewise  linear,  and  (6)  piecewise 
s i nus o i da 1 / p i e cew i s e sinusoidal. 

Whenever  the  wire  structure  is  of  sufficiently  large 
magnitude,  in  terms  of  numbers  of  wavelengths  of  wire,  one  may 
not  be  able  to  require  the  increment  A to  be  so  small  relative 
to  wavelength  to  be  assured  that  piecewise  sinusoidal  functions 
(for  basis  or  testing  set)  yield  results  which  are  extremely 
close  to  those  obtained  by  means  of  piecewise  linear  functions. 

For  A/A  ^1/10  the  difference  should  be  less  than  6 or  7 percent. 

The  solution  procedure  outlined  above  in  which  piecewise 
linear  functions,  or  triangles,  are  utilized  as  elements  of 
both  the  basis  and  testing  sets  is  an  example  of  a Galerkin 
method,  i.e.,  a method  in  which  testing  functions  are  the  same 
as  basis  functions.  Since  (2)  and  (6)  above  are  approximately 
equivalent,  the  procedure  employing  triangles  for  expansion  and 
testing  is  quite  close  to  Galerkin's  method  with  piecewise 
sinusoidal  functions. 
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2.8  ANALYSIS  OF  STRAIGHT-WIRE  ANTENNA 


The  solution  methods  discussed  for  the  scatterer  in 
Section  2.7  are  directly  extendable  to  the  case  of  a straight- 
wire  antenna,  driven  by  a voltage  source  of  V volts  and 


located  at  some  point  along  the  length  of  the  wire  (Fig.  2.2) 

If  one  adopts  the  standard  slice  generator  source  model  [4], 

tlien  E^  of  (2.1)  becomes 
z 


E (z)  = V6(z-z  ) 

z g 

where  6(z)  is  the  familiar  delta  function.  Without  the 


(2.32) 


proximations  of  Section  2,4,  V of  (2.21)  becomes  exactly 

m 


— V I 6(z-z  )A^’^(z)dz 
w I g m ' " 


-j— VA^’®(z  ) 

OJ  m g 


(2.33) 


for  either  piecewise  linear  or  piecewise  sinusoidal  testing 
and,  for  either  testing  set,  reduces  to 


V = -j V , (z  =z  ) 

m 0)  g m 


(2.34) 


when  z^  and  z^  coincide.  For  multiple  excitation  (2.33) 
becomes 


(2.35) 


The  remaininB  terms  in  (2.21),  S and  I , are  handled 

ran  n 

in  exactly  the  same  way  that  they  are  for  a scatterer.  Hence, 
with  the  appropriate  definition  of  in  (2.33)--(2.35)  and 
with  for  pulse  or  triangle  representation  of  current  dis- 

cussed in  Section  2.7,  one  may  calculate  current  on  an  arbitrarily 
excited,  s t r ai gh t - w i r e antenna. 

Sample  results  for  a center-driven,  hal f-wave length 
antenna  are  depicted  in  Fig.  2.8.  Current,  calculated  by 
means  of  piecewise  linear  expansion  and  testing,  is  seen  in 
Fig.  2.8,  where  for  comparison  is  also  displayed  current  cal- 
culated from  a pulse  expansion  with  collocation  numerical 
solution  to  Hallen's  equation  [5]. 
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Hallen  ( 36x36 -Pulse/Collocation) 

xxxx  Eq.  2.1  (15x15 -Triangle/Triangle) 


Fig.  2.8.  Current  on  Center-Driven  Dipole 
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CHAPTER  3 


ANALYSIS  OF  COUPLED  WIRES 


3.1  INTRODUCTION 


With  tlie  foundation  being  established  in  Chapter  2,  the 
next  step  in  the  procedure  for  developing  a general  analytical 
technique  for  characterizing  wire  structures  is  the  treatment 
of  a particular  wire  in  the  presence  of  others.  The  electro- 
magnetic influence  on  one  wire  in  an  ensemble  due  to  all  others 
plus  itself  must  be  fully  accounted  for  in  any  analysis.  Such 
coupling,  mutual  and  self,  must  be  taken  into  consideration  for 
each  wire  in  an  ensemble. 

In  this  Chapter  .are  presented  the  analysis  of  wire  struc- 
tures and  the  formulation  of  pertinent  integral  equations. 
Numerical  solution  procedures  are  discussed  as  are  s y mme  try 
conditions.  To  support  the  theory,  results  of  sample  calcula- 
tions are  included  and  compared  with  data  obtained  elsewhere. 

In  Fig.  3.1  are  shown  several  wires,  one  designated  the 
p^^'  wire  and  another  designated  the  q^*^.  We  focus  major  atten- 
tion upon  the  p^^  wire  and  investigate  how  one  may  calculate 

the  electric  field  on  this  wire  due  to  the  current  i on  itself 

P 

plus  that  due  to  i , the  current  on  the  q wire.  The  contri- 

q . 

bution  on  the  p*”'^  wire  due  to  the  presence  of  wires  other  than 
the  q*"^  may  be  calculated  in  an  identical  way,  and  then  the  total 


contribution  from  wires  otlier  than  the  p is  simply  a summation  over 


V • 


all  wires,  q = 1 , 


2,...,  N(q^p).  Written  symbolically,  the 


.•lectrir  tielcl  directed  along  the  axis  ol  the  p wire  and 
evaluated  on  its  surface  is 


E^(0 


= E (C)  + u 


E (?) 


(3.1) 


q = l 

(q?^p) 


where  is  the  total  axially  directed  field  on  the  surface  of 

the  p^^*  wire  produced  by  all  charges  and  currents  on  the  struc- 
ture, E is  that  due  to  the  current  i , and  E is  that  due  to 

p q 

the  wire.  In  Eq.  (3.1)  u^  is  the  unit  vector  along  the 

p*”^  wire  (and  has  the  sense  of  the  current  i ),  and  ? is  an 

P 

independent  variable  denoting  axial  displacement  along  u^  of 
the  p*"^  wire  in  a local  coordinate  system  with  origin  at  the 
center  of  the  p*"'^  wire.  Fig.  3.2  depicts  the  two  wires  and 
serves  to  define  geometric  quantities  of  interest.* 

At  a general  point  in  space  (x,y,z)  one  may  write  E 
in  the  following  form: 


E (r)  = 


-jo:Ap(r)  - $p(r) 


(3.2) 


where  A^  and  are  the  vector  and  scalar  potentials,  respectively, 

t h 

calculated  from  the  sources  on  the  p wire.  Subject  to  the 
thin-wire  approximations,  these  potentials  may  be  written 


*See  Appendix  III. 


^ Vt  ' --  - r 


-jkl 


= Air 


i (A’) 
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— d^' 


r'=-L  /2 
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f J . 3ai 


L /2 
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$ (r) 
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^ f 

rire  J 


C'=-L  /2 
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P,  (C) 
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-jk I r-r  ' 

. P 


r-r’  [ 
P 


+ t (r) 

P P 


(3.3b) 


where  the  general  point  (x,y,z)  is  located  by 


r=xu  +yu  +ZU 
X y z 


(3.  A) 


and  a source  point  on  the  p wire  by 


r ' = r + Z'  u 

P P P 


(3.5) 


“ C t ll 

In  (3.5),  r^  is  a vector  wliich  locates*  the  center  of  the  p 

wire  so  (3.5)  itself  is  simply  the  equation  of  a straight  line 

in  space  along  the  axis  of  the  p^^  wire.  The  linear  charge 

density  on  the  p*"^  wire  is  denoted  p^^^  and  is,  of  course,  related 

P 

to  the  p*"*^  wire  current  by  the  continuity  equation. 


*See  Appendix  III 


dC 


i (O  + jojp.  U)  = 0 


( 'J  . 6 > 


r + r ~ 

In  (3.3)  the  t.wo  terms,  and  , represent  contributions  to 

the  scalar  potential  from  the  rings  of  discrete  charge  located, 
respectively,  at  the  upper  and  lower  ends  of  tiie  p*"''  wire. 

These  charges  are  present  whenever  the  current  is  not  zero  at 
tl>e  ends  as  is  necessary,  if  the  analysis  is  to  be  extended  to 
structures  involving  wires  which  join  at  their  end-points. 


th 


Similarly,  tlie  electric  field  E^(r)  due  to  sources  on  the 
wire  may  be  written 


E (r)  = -jcjA  (r)u  - grad  <I>  (3.7) 

q q q q 

wliere  ®'^d  are  the  vector  and  scalar  potentials,  respec- 
tively, calculated  from  the  sources  on  the  wire.  A and 

q 

can  be  expressed  in  terms  of  potential  integrals  like  those 
of  (3.3)  with  p replaced  by  q.  The  unit  vector  u^  is  directed 
along  the  q j.  wire  and  has  the  sense  of  the  defined  current  on 
this  wire. 

Since  tiie  contributions  on  the  p^^  wire  due  to  currents  on 
tlie  other  wires  are  all  of  the  form  of  (3.7),  only  two  wires 
are  treated  in  general  below  (Fig. 3. 3)  so  that  the  discussion  is 
not  rendered  overly  complex.  The  extension  to  more  than  two 
wires  is  quite  direct  and  can  be  done  readily  after  the  present 
two-wire  case  foundation  has  been  established. 


4 5 


llif  .ixially-dirfcted  electric  field  on  the  surface  of 

t h 1 1 1 i 

Lin.'  p wire  siiown  in  Fi^-  '3.2,  due  to  sources  on  botti  Liie  p 

. t h . 

and  q wires,  is 


h.  = -juA  - V-  $ 

C P 3 C P 


i u)A  ( u • u ) - 7^  <5> 

q q P q 


(3.8) 


Retaining  but  employing  the  Lorentz  gauge  explicity  to 

eliminate  $ one  may  convert  (3.8)  to  an  alternate  form: 

P 


^ + kM  A (o 
0^2  J P 


2 

-f-  k^A  (0(u-u  ) - 3^  I7  $ (O 
q q p U dC  q 


(3.9) 


Actually,  the  Lorentz  gauge  is  invoked  implicity  relative  to 

for  otherwise  the  potential  integral  representation  of  the 

scalar  potential  would  be  inappropriate;  however,  in  (3.9)  the 

relationship  between  and  A^  is  not  invoked  directly.  In 

view  of  the  he te rogeneous  application  above  of  the  relationship 

between  vector  and  scalar  potential  under  the  Lorentz  gauge, 

one  must  exercise  care  to  avoid  violation  of  the  continuity 

t h 

equation,  when  he  calculates  from  the  charge  on  the  q 

wire.  When  this  wire  joins  others,  as  is  anticipated  by  the 


expression  (3.3b),  special  heed  must  be  given  this  caution. 


W' 


and  one  must  recall  that  explicit  application  of  the  Lorentz 
gauge  ensures  a complete  electric  system  [IJ. 


3.2  TESTING 

As  is  done  above  for  the  isolated  wire  equations,  we  seek 

2,  s 

a set  of  linear  equations  by  testing  with  and  of  (2.6) 

pursuant  to  the  development  of  a numerical  procedure  for  per- 
forming various  calculations  on  wire  structures.  As  before  we 
distinguish  between  interior  testing  and  end  testing  in  order 
that  tlie  boundary  conditions  on  wire  ends  may  receive  desired 
attention  and  may  be  efficiently  incorporated  in  the  numerical 
solution  technique. 

Testing  of  (3.9)  with  a testing  function  T yields 


wtiicli,  in  view  of  the  definition  of  the  scalar  product  (2.4b) 
and  the  subdomain  nature  of  the  testing  functions  to  be  con- 
sidered here,  can  be  written 
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/•  > /s  /V  /* 


A (i;;  T (c)d'- 
q m 


(C)  T (OdC 

d C d ® 


E_  (C)T^(Od; 
C m 


(3.11) 


As  before,  in  the  case  of  end  testing,  the  limits  of  the  above 
integrals  must  be  modified;  in  particular,  for  testing  over  the 
lower  interval,  (C  ^ M+1^’  lower  and  upper  limits  of 

each  integral  are  and  ^ > respectively,  whereas,  for 

testing  over  the  upper  interval,  CE  respective 

lower  and  upper  limits  are  2 course,  one  recog- 


nizes that  C.,,  = IL  1 2 and  that 

1 M p ± 


(M-1  ) 


/ 2 - A ),  where 
P P 


is  the  subdomain  interval  length  defined  in  Fig.  3.3  . 

Interior  Testing  Making  use  of  the  equations  in  Appendix  I, 
one  can  perform  the  operations  indicated  in  (3.11)  to  achieve 
the  equations  resulting  from  testing  with  either  piecewise  linear 

or  piecewise  sinusoids  over  the  interior  intervals.  For  piece- 

2, 

wise  linear  testing  (T^  = A^)  one  obtains  Eq.  (3.12)  and  for 

piecewise  sinusoidal  testing  (T  = A ) he  obtains  Eq . (3.13), 

mm  ^ 

both  of  which  are  found  on  the  following  page.  Eqs.  (3.12)  and 
(3.13)  are  exact  in  the  sense  that  the  testing  procedure  has 
Introduced  no  approximations. 


s e Line 


I 


Kiid  In  end  testing  we  consider  bolh  Kq  . ( i . 8 ) and 


s y- 

[•:(|.  (i.9)  and  test  witli  cacti  of  the  functions  A , and  A ^ . 01 


irse,  these  equations  are  equivalent  and  either  may  lie  used 


in  any  wire  analysis,  or  they  both  may  be  employed  in  a mixed 


fashion.  Testing  Eq . (3.8)  by  means  of  over  the  upper  and 


lower  end  Intervals  yields  (3.14a)  and  (3.14b),  respectively, 


and  end  testing  (3.8)  with  yields  (3.15a)  and  (3.15b) 


Similarly,  testing  Eq.  (3.9)  with  A^j^  over  the  upper  and  lower 


end  intervals  yields  (3.16a)  and  (3.16b),  respectively,  while 


end  testing  of  (3.9)  with  A^j_j  gives  rise  to  (3.17a)  and  (3.17b). 


One  should  note  tliat,  in  all  end  testing  cases  mentioned  here. 


the  scalar  product  defined  in  (2.4b)  is  employed.  The  limits 


on  the  r i gh t -h and-s ide  integral  in  (2.4b)  applied  to  the  present 


case  would  be  -L^/2  and  +L^/2,  which  reduce  end  testing  to 


integration  over  only  one-half  of  the  domains  of  the  testing 


£ S 

functions,  or  this  is  the  reason  that  end  testing  is 


treated  apart  from  interior  testing. 


(?)dC:  (3.14b) 


e a 


3.3  APPROXIMATIONS 


As  was  tile  case  in  the  analysis  of  Llie  Isolated  wire  of 
Chapter  2,  there  are  good  approximations  whicli  can  be  used  to 
advantage  to  significantly  simplify  the  present  expressions. 
The  integral  approximations. 


and 


/ 


-M+1 

f(C)A^’"(Od?  = i Apf(C_„) 


/ 


f(OA^’"(OdC  = i Apf(^^) 


C = C 


M-1 


(3. I8a) 


(3.  28b) 


are  good  wlienever  f(^)  is  smooth  over  the  lower  end  interval  in 
the  case  of  (3.18a)  and  over  the  upper  end  interval  in  the  case 
of  (3.18b).  These  approximations,  (2.15)  and  (3.18),  are  applied 
in  (3.12)  - (3.17)  to  integrals  containing  the  electric  field 

and  vector  potential  terms,  A^  and  A^,  and  one  thereby  arrives 
at  the  following  corresponding  approximate  equations  (see  fol- 
low ing  pages). 
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Piecewise  Linear  End  Testing  with  Approximation  (Case  I - Te s ting  of  ( 3.8)) 


(3.22b) 


Other  approximations  may  be  made  in  the  present  analysis 
but  they  are  ol  a type  wliose  error  depends  on  parameters  in 
the  caieulations  b e i n p performed.  For  example,  if  the  deriva- 
tive ot  a lunction  f ( C ) is  sufficiently  smooth  over  the 

range  of  integration  indicated  below,  one  may  employ  tiie 
approximations , 


^in+1 


(3.23a) 


, 2.  ,s 


(C)dC 


— ^ — f (C  ) 

2 3;; 


and 


I 


M+1 

3 


A^’"(?)dC 


-JS. 

2 3c 


-M 


f (C  m) 


(3.23b) 


(3.23c) 


Approximations  (3.23)  may  be  applied  to  terms  like  tlie 

third  in(3.1i)  Involving  <I>^  to  arrive  at  simpler  expressions 

than  those  appearing  in  (3.12)  - (3.22).  Such  approximations 

are  of  acceptable  accuracy  only  wlien  is  sufficiently 

smooth  over  the  range  of  integration  of  the  pertinent  expres- 

Sion  (3.23);  ?—  $ is  smooth  enough  for  (3.23)  to  apply  in 

3 C q 
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•S  i lu.it  i Uiis  where  ^ is  the  scalar  potential  due  to  an  clement 
ot  current  remote  (In  terms  of  wavelength)  from  the  interval 
over  which  one  of  the  approximations  (3.23)  is  to  be  employed. 


3.4  USE  OF  PIECEWISE  LINEAR  TESTING  IN 
NUMERICAL  SOLUTION  PROCEDURE  FOR  DETERMINATION  OF 
CURRENT  ON  SCATTERER  COMPRISING  TWO  COUPLED  WIRES 


In  this  section  is  outlined  the  theory  for  calculating 
currents  on  a two-wire  scatterer  illuminated  by  a specified 
incident  field.  So  that  explanations  are  not  overly  complicated, 
attention  here  is  limited  to  two  wires  but  the  technique  is 
general  and  can  be  applied,  with  little  or  no  modification,  to 
structures  of  many  wires.  Pulses  are  used  for  representing  the 
current  even  though  other  basis  sets  can  be  employed,  if 
desirous . 


As  is  done  in  Chapter  2 for  the  straight  wire,  we  first 
approximate  the  current  on  each  of  the  two  wires  by 


P 

!■ 


i„(z)  = ^ij’i  fz) 
n = -N 


(3.24a) 


and 


i^(z)  = 


P 

2 


I^ir,  (z) 

n n 


(3.24b) 


n = -N 


where  (2N  +1)  and  (2N  +1)  are  the  numbers  of  terms  in  the 

P q 

approximation  for  i^  and  i^,  respectively.  Next  tlie  approxima- 


t i on  (3.24a)  is  substituted  into  Kquations  (3.3)  and  (3.24b), 

into  a similar  pair  for  A and  . With  tlie  currents  in  tiie 

q q 

potential  integrals  replaced  by  (3.24),  A , '1'  , A , and  'll 

p p q q 

are  employed  in  the  selected  tested  equations--either  (i.ly) 
and  (3.21)  or  (3. 20)  and  (3.22)--  which  sequence  of  steps  esta- 
blishes a system  of  linear  equations  having  many  of  the  features 
of  (2.21).  In  this  study  of  two  wires,  and  subsequent  treat- 
ments of  multiple  wires,  the  wire  currents  are  approximated  by 
pulses,  and  piecewise  linear  functions  are  used  for  testing 
(Equations  (3.ly)  and  (3.21)). 

In  preparation  for  writing  (3.19)  and  (3.2l)  as  functions 
of  the  and  it  is  desirable  to  develop  expressions 

for  the  vector  and  scalar  potential  terms  appearing  in  these 

two  equations.  The  vector  potential  A on  the  wire  due  to 

P 

current  i^  on  that  wire,  differs  from  the  corresponding  quantity 
in  Chapter  2 only  in  that  contributions  from  the  end  half-pulses 
must  be  included.  One  can  readily  show  that  at  a point  ^ on 
the  p wire  (in  the  local  coordinates  of  the  p*"^ 
origin  at  its  center) 


wire  with 
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(3.25) 


where  the  first  and  last  integrals  are  contributions  from  the 
lower-and  upper-end  half-pulses,  respectively,  and  where 


^ =mA,rn=0,±l,  ±2,... 

m p 


(3.26) 


> 


i 

V 


The  scalar  potential  on  the  p wire  due  to  the  charge  on  this 

wire  can  be  calculated  with  the  aid  of  Fig.  3.4  and  Appendix  II. 

t h 

As  suggested  in  Fig.  3.4,  associated  with  the  n current  pulse 

of  strength  is  a ring  charge  at  ^^-A^/2  having  total  charge 

-I^/ioj  and  another  at  C +A  /2  having  total  charge  +I*^/jcj.  The 
n n p n 

total  quantities  of  charge  of  the  ring  charges  deposited  by  the 
lower-end  pulse  are  +I^j^  at  -L  /2  and  /jw  at  “L  /2  + A /2 

p p - p p p 

and  those  deposited  by  the  upper-end  pulse  are  /jw  at 

P 

L /2  - A /2  and  if,  /iw  a*  L /2.  Therefore,  the  scalar  potential 
P P Np  p 
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• U .1  Joiaiiou  f,  on  liio  p wire  due  to  i L .s  own  charges  is 


K ( C+A  [ N - ] ) - K ( ;+A  N 

1>  P PI 


N -1 
P 


n=-(N  -1) 
P 


K(C-Ap[n-  >-i  ] ) 


+ L^C-ApNp)  - K(C-Ap[Np-  J,]) 


(3.27  ) 


Ktiuations  (3.25)  and  (3.27)  provide  a means  of  displaying  tiie 
dependence  of  quantities  bearing  subscript  p upon  in  ( 3 . l9  j 

and  (3.2  ).  Quantities  bearing  subscript  q in  these  two  equations 
are  next  written  with  explicit  dependence  upon  I^.  Involved  in 
eacii  such  quantity  is  the  distance  from  an  appropriate  source 
point  on  the  q*"^  wire  to  a field  point  on  the  surface  of  the  p*"*^ 
wire.  This  distance  D is  derived  in  Appendix  III  where  geo- 

pq  6 

metric  quantities  are  defined  and  discussed  in  detail: 


D (s 

pq 


4.S')  - [aj  + Kf'a-Up)  - 


■vil'' 


(3.28) 


The  vector  potential  at  a point  on  the  p wire  due 

to  the  current  (3.24b)  on  the  q^*^  wire  is 
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-(N  - 3s)A 

q q 


A (C  ) = 7^/ 
q m 4ti  ^ -N 


/ 


G (A  , A ’ ) d C ' 
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C'=-N  A 

q q 


N -1 

q 


(n+  ) A 


2 ■:  / 


n--(N  -1)  C'=(n-  u )A 


N A 

q q 


+ 1 


I 


(3.29) 


wh  ere 


-jkDpq(;.C') 

^ ~ Dp/T;T') 


(3. 30) 


The  scalar  potential  is  formulated  in  a manner  similar 
to  that  leading  to  (3.27)  but  with  the  kernel  (3.30)  playing  the 
role  ofthe  usual  s t r a i gh t - wi r e kernel.  Using  (II-7)  of  Appendix  II, 
one  can  show  that 


i tr  ) = k ...  ) \c  (C  , [-N  + I5] A ) - G (;  ,-N  A ) 

^ i^incueN  -N^L  P‘5  ‘5  ^ ‘5  P‘5  ‘5  q. 


q 

+ ^ I*’  fc  (C  ,[n+  !i]A  ) - G (C  ,[n-  ] A ) 

/ n 1_  pq  q pq  q J 


n=-(N  -1) 

q 


[Sq^^ 


(3.31) 


Linear  System  of  Equations  Witli  A of  (3.25)  and  <I>  of  (3.27) 
written  in  terms  of  the  unknown  coefficients  and,  also, 

with  A of  (3.29)  and  $ of  (3.31)  written  in  terms  of  the 

q q 

I^'s,  one  substitutes  these  vector  and  scalar  potential  quan- 
tities into  the  tested  equations  (3.19)  and  (3.21).  As  a result 


with  M=Np , one  has  a system  of  linear  algebraic  equations  like 

(2.21)  , 


I^S  + 
n mn 


I'^G  = V 


m = 0 , ±1  , ±2 , + , . . . ,±N  , 

P 


(3.32) 


which  is  a system  of  2N  +1  equations  in  2(N  +N  +1)  unknowns. 

P P q 

The  coefficients  S are  called  "self  terms"  and  renre- 

mn 

sent  the  entries  in  (3.19)  and  (3.21),  subject  to  replacement 
of  quantities  by  their  equivalents  in  (3.25)  and  (3.27),  which 


«NII  . 


V 


,1  r *•  f (>(■  1 t i c i *Mi  t s () ! the  i’  u r r e n L h 1 ^ 

n 


S 


n 


/]  n (1  S 


a r V 


11  ‘ ' - N , n ' , 11 

P P 

c 111.-  1 ! i c i e 11 1 * of  I n = - N , . . . , 0 , . . . , N , in  ( J . 2 1 a > and  ( d . 2 J b ) , 

11  p p 

respectively.  S is  the  coefficient  of  in  (3.19),  m=-(N  -1), 
' mil  n p ’ 

...,0,...,N  -1.  The  coefficients  C in  (3.32)  are  called  "cou- 
p mn 

pling  terms"  since  they  are  associated  with  contributions  from 

the  current  on  the  wire  to  vector  and  scalar  potential 

quantities  evaluated  on  the  p*"^'  wire.  These  coefficients  may 

be  readily  identified,  if  one  substitutes  (3.29)  and  (3.31) 

into  (3.19)  and  (3.21).  C and  C.,  , n = -N  ,...,0,...,N  , 

-N,n  N,n’  q’  q 

P P 

are  coefficients  of  in  (3.21a)  and  (3.21b),  respectively, 

while  C ,m  = -(N  -1),...,0,...,N  -1,  is  the  coefficient  of 

m , n P P n 

in  (3.19).  In  other  words,  the  S 's  account  for  contributions 

mn 

to  vector  and  scalar  potential  at  locations  along  the  p ^ wire 

from  the  individual  elements  of  the  basis  set  representing  i , 

whereas  the  C 's  account  for  contributions  on  the  wire 

mn  ^ 

due  to  elements  of  the  basis  set  by  means  of  which  i is 

q 

approximated . 

The  inhomogeneous  or  driving  terms  V 's,  m=-N  ,...,0,...N  , 

m p p 


(3 . 33a) 


*Commas  are  used  to  separate  subscripts  when  their  being  in 
juxtaposition  may  lead  to  misinterpretation;  S and  S are 

..  u 1 1 j .-1  1 

to  be  looked  upon  as  the  same. 


(3.33b) 


V 

ni 


-A  E (C  )-u 
P m p 


m=- (N  -1)  , . . . ,N  -1 
P P 


c'i  n d 


(3. 33c) 


which  are  the  right-hand  side  terms  in  (3.19)  and  (3.21)  with 

E (C  ) replaced  by  -E^(C  )*u  so  that  boundary  conditions  or 

^ m ' m p 

the  surface  of  the  wire  are  enforced.  In  the  case  of  an  antenna 
problem,  the  modification  explained  in  Section  2.8  is  utilized. 

The  system  of  equations  (3.32)  is  general  in  the  senses 
chat  the  p*"^’  wire  and  the  q''*'  wire  may  be  any  pair  in  a com- 
plex wire  structure  and  that  their  relative  orientation  is 
completely  arbitrary  (Appendix  III).  To  facilitate  subsequent 
discussions  we  write  (3.32)  as  a matrix  equation, 

[Sp]  [iP]  + [cj]  [l'’]  = IVp]  , (3.34) 

and  define  the  individual  matrices  below. 
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(3. 35e) 

The  subscript  p on  the  matrices,  [S  ],  [C*^],  and  [V  ],  implies 

P P P 

that  each  of  these  is  associated  with  evaluation  of  quantities 

on  the  wire.  The  superscript  p on  [I^]  identifies  this 

matrix  as  the  column  vector  of  current  coefficients,  1^'s, 

n 

in  the  expansion  of  i^  on  the  p*"*^  wire;  similarly,  [l'']  is  the 
vector  of  coefficients  in  the  approximation  of  i^.  The  super- 
script q on  [Cp]  denotes  that  this  is  a matrix  of  coupling 
terms  indicating  a contribution  from  the  q*'*^  wire  current. 

Equation  (3.32),  or  its  matrix  equivalent  (3.3A),  serves 
to  enforce  the  boundary  condition  that  the  axial  electric  field  on 
the  p^^'  wire  be  zero.  Of  course,  this  boundary  condition 
must  be  satisfied  on  all  wires  in  a given  structure,  and  (3.32) 
and  (3.34)  are  totally  general  so  that,  by  a mere  change  of 
subscripts  and  superscripts,  one  has  available  the  needed 
equations  to  enforce  this  boundary  condition  on  other  wires. 

For  example,  if  two  wires  which  are  located  and  oriented  in 
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1 


i 


cyan  , j^j^>  ^ ^ m^^pm  M - iC^  w4 


■^ivfiilfi 


sonu-  specified  way  (Appendix  III),  are  illuininaLed  by  an 
incident  lie  Id,  the  entire  pro  1)  Jem  can  be  cliaracterized  tjy 
tile  pair  of  matrix  equations: 


[Sp]  [iP]  + [Cj]  [I^]  = [Vp]  , (3.36a) 

[cP]  [jP]  + [S^j  [I'l]  = [V^]  . (3.36b) 

Tiiese  suggest  that  the  wires  are  designated  and  and 

(3.36a)  is  recognized  to  be  (3.3A);  (3.36b)  is  the  companion 

equation  whose  matrix  elements  are  determined  from  the  same 
procedure  used  to  calculate  those  in(3.36a)  (with  indices 
interchanged).  If  there  were  a third  wire  designated  the 
w^*^  wire  in  the  problem  under  consideration,  the  appropriate 
equations  would  be 


[IP] 

-f- 

[Cj] 

+ 

[l'"]  = 

[cjj 

[IP] 

+ 

[s'^] 

[l'’] 

+ 

[I”]  = 

[V  ] 

q 

[IP] 

+ 

t^wl 

[s'^l 

[l“]  = 

or,  as  a single  matrix  equation. 


r 

- 

- 

[S  ] [c^]  [c“] 

[IP] 

[V  ] 

p p p 

p 

[cP]  [SI  [c“] 

[l'’] 

_ 

q q q 

q 

[cP]  [C,^]  [S^] 

[I'"] 

WWW 

W 

. 

- 
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(3.37) 
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Wires  Not  .)  o 1 ae  d If  tlie  wires  in  a sLructure  under  stuu. 

.ire  iioL  in  electrical  contact,  the  current  at  the  free  wire 
ends  must  he  zero.  Therefore,  t li  e coefficients  of  tlie  end 
current  pulses,  i.e.,  1^^^  and  1^  tor  the  p*"*'  wire,  are  made  equal 
to  zero  (actually  discarded  from  the  set  of  unknown  coefficients) 
Discarding  these  upper-  and  lower-end  coefficients  reduces 
the  number  of  unknowns  so  one  in  turn  discards  the  equations 
in  the  system  of  linear  equations  (3.37)  which  correspond  to 
end  testing. 

If,  for  example,  the  wires  in  the  two-wire  problem  char- 
acterized by  (3.36)  are  not  joined  electrically,  one  would 

set  , ij!  , and  1^,  equal  to  zero.  Also,  he  would 

- .N  N , - N N 

p p p q 

delete  those  equations  resulting  from  end  testing.  That  is  , 


the  following  matrix  elements  would  he  deleted:  S 


N ,n’ 
P 


and 

S^,n’ 

n = -N 

p 

and 

,n’ 

P 

n = -.V 

q 

-N  ,n  ’ 

q 


-N  ’ N 
P P 


-N^,  and  . After  these  deletions,  the  resulting  linear 

q 

system  would  contain  2 (N  +N  -1)  equations  and,  of  course, 

P q 

the  same  number  of  unknowns. 

Two  Wires  Joined  As  an  example  of  a simple  wire  structure 
having  a junction,  we  consider  the  two  wires  illustrated  in 
Fig.  3.5  where  one  sees  that  the  upper  end  of  the  p''^’  wire 
joins  the  lower  end  of  the  q^  wire.  Of  course  at  the  free 
ends,  lower  end  of  the  p*^^  wire  and  upper  end  of  the  q^^. 


ihe  lia  If -pulse  currents  are  set  equal  to  zero,  and  the  matrix 


equation  is  modified,  as  described  above,  to  be  consistent 

with  i*^j^  since  the  current  must  be  continuous 

* P <1 

at  the  iunction,  one  requires  tliat  Ij^  be  equal  to  1^,  , which, 

l>  ''I 

ol  course,  reduces  (lie  nuinlu'.i,  ol  unknowns  In  Ihe  s y s I e in  nl 

equations  by  one,  thereby  requiring  further  modification  <jf 

the  matrix  equations  (3.3b). 

With  if.  (junction  current),  the  columns  below 

N -N  J 
p q 

both  operate  on  Ij  (recall  C are  zero 

P P q q ’ 

deleted  from  matrix--to  account  for  =1?,  =0  at  the  free 

-N  N 

p q 

ends  ) : 


” 

^-N  +1 ,N 

fl.-N 

P P 

p q 

S 

c „ 

N , N 

N ,-N 

P P 

(3.38a)  , 

p q 

(3.38b) 


Hence  Llie  above  two  columns  are  added  to  reduce  the  number  of 
columns  by  one.  Next  the  number  of  rows  must  be  reduced  by 
one  too.  Reduction  of  the  number  of  rows  by  one  is  achieved 
by  adding  the  and  rows  on  both  sides  of  the  equation, 

The  summed  row  on  the  left-hand  side  is 


,-N 

P 


+ 1 ‘'-N  ,-N  +1 

P q P 


^N  ,N  ,-N  ■*"  ^-N  ,N  ^-N  ,-N 

pp  pq  qp  q q 


^N  ,N  -1  ^-N  ,N  -1 

p q q q 


(3.39) 


while  that  on  the  right-hand  side  is  simply 


V + V 

N -N 

L p q . 


(3.40) 


The  center  entry  in  (3.39)  comprises  four  terms  because  these 
terms  are  members  of  the  summed  column  vectors  (3.38)  as  well 
as  of  the  summed  row  vectors. 

As  an  alternative  to  the  first  step  in  the  modification 
described  above,  one  may  elect  to  append  to  the  system  (3.36) 
an  additional  equation  enforcing  continuity  of  current  at  the 
junction,  1^^  » rather  than  to  add  the  column  vectors 

^P  ~ q 

(3.38).  Of  course,  the  second  step  of  adding  rows  must  be 
done  regardless  of  how  the  current  continuity  is  handled. 


r 


With  the  matrix  modifications  discussed  above,  the  re- 
sulting linear  system  of  equations  possesses  the  same  number 
of  equations  and  unknowns,  and,  therefore,  one  may  employ 
standard  methods  to  determine  the  unknown  current  coefficients, 

I^'s  and  I'^'s. 
n n 

Sample  Results  As  examples  of  the  utility  of  the  coupled 
thin-wire  theory  and  numerical  techniques  outlined  in  this 
Chapter,  calculated  currents  on  bent-wir^  (joined)  scatterers 
are  presented  in  Figs.  3.6,  3.7,  and  3.8.  Wire  configuration 

and  dimensions,  as  well  as  excitation,  are  indicated  in  the 
figures,  and, where  available,  results  from  other  calculations 
are  provided  for  comparison. 

Extensive  data  calculated  by  means  of  the  theory  given 
in  Chapters  2 and  3 can  be  found  in  Chapter  4. 


I 

1 

i 


Dependence  of  Current  on  Bent-Wire  Scatterer  on  Bend  Angle 


CHAPTER  4 

APPLICATION  TO  WIRE  BICONICAL  ANTENNA  STRUCTURES 


4.1  MATRIX  EQUATIONS  AND  SYMMETRY 


The  coupled  integral  equations  formulated  in  Chapter  '3, 
along  with  either  piecewise  linear  testing  or  piecewise  sinu- 
soidal testing,  leads  one  to  generalized  equations  (3.12)  to 
(3.15)  for  analysis  of  wire  structures.  Since  the  pulse  basis 
set  is  easy  to  incorporate  and  is  as  powerful  as  any  other  set, 
only  tlie  pulses  are  used  in  the  following  application  to  bi- 
conical  antennas.  Sample  results,  viz  .,  the  current  distribution, 
input  impedance  at  the  feed  point,  and  the  far  electric  field 
radiation  pattern,  are  presented  for  a wire  biconical  antenna 
in  free  space  or  above  a ground  plane  in  various  orientations. 

From  the  developments  in  Chapter  3,  a self  matrix  [S^]  and  a 

P 

mutual  coupling  matrix  [ C*^  ] are  established  for  two  arbitrarilv 

P ^ 

oriented  wires,  p*"^  and  either  joined  at  one  end  or  not. 

For  the  case  of  a biconical  antenna  above  a ground  plane. 

Fig.  4.1,  or  a biconical  antenna  in  free  space.  Fig.  4.2,  the 
following  partitioned  matrix  equation  of  eight  by  eight  is 
obtained , 
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[sj]  [c^l  [C^]  [cj] 


[S2]  [C^]  [cj] 


[sP]  [C^] 

p p 


[c'l  (C^l  [C^l  lc‘l 


(4.1) 


12  3 8 

where  [I  ],  (I  ],  [I  ] are  partitioned  column  vectors 

representing  the  unknown  current  coefficients  on  the  first, 


second,  ..., eighth  wires,  and  [V^],  [V2],  are 


the  corre- 


sponding excitation  column  vectors  defined  by  (3.35). 

The  self  and  the  mutual  coupling  partitioned  matrices, 

[Sp]  and  [Cp],  result  from  piecewise  linear  testing  on  the  p*”^ 
wire  and  account  for  contributions  from  the  currents  on  p*"^ 
and  wires,  res  pe  c t i ve  ly  . Due  to  the  various  symmetries* 

in  the  orientation  of  the  wires  and  in  the  location  of  the  feed 
point  at  the  center  of  the  structure,  [I^]  = [I^],  = [I^], 

[I^]  = [I^]  and  [I^]  = [I^].  The  bottom  half  of  the  matrix 

can  be  generated  using  top-half  matrix  elements  in  the  reverse 
order.  Hence,  deleting  the  bottom  half  of  the  partitioned 
matrix  (4.1),  one  needs  to  generate  only  the  top  half  of  (4.1): 


*These  column  vectors  are  equal  subject  to  appropriate  ordering 
of  the  elements  in  the  columns. 


JDK,-? 
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[sjl 

[cj]  [cj]  [c\] 

• [C®] 

[1^1 

[Vil 

[C^] 

[s2] 

[C^] 

• [C®] 

[I^] 

[V2] 

[C3] 

[C3] 

[S3] 

• [C®] 

[I^] 

[V3] 

_[ci] 

[C^] 

[C^]  ■ 

• [c®] 

[I^] 

(4.2) 

Due  to  the  symmetry  in  the  currents  as  stated  above,  the  par- 
titioned matrix  equation  (4.2)  can  be  folded  and  written  as 


[S J+c®] 

[C^i+cj] 

[C3+C®] 

4 5 

[Ci+Ci] 

[I^] 

[ViJ 

[cj+c®] 

2 7 

IS2+C2 ] 

[C2+C2] 

[C2+C2] 

[I^J 

[V2] 

[C3+C3] 

2 7 

[C3+C3] 

u^] 

[V3] 

[cj+c®] 

[C^+C^] 

ic]+c‘l 

IsJ+Cj) 

u^] 

_ 

[V,] 

(4.3) 


The  partitioned  matrix  equation  (4.3)  is  a square  matrix  of 

four  by  four  and  can  be  solved  for  the  unknown  current  vectors 
12  3 4 

[I  ]»  [I  ]>  [I  ],  and  [I  ].  This  compact  matrix  form  is  equi- 
valent to  successive  modifications  of  the  kernel  functions  in 
the  integral  equations  (3.12)  and  (3.14).  Such  are  possible 
because  of  the  symmetry  in  the  d 1 s t r 1 b u t i on o f the  currents 

with  respect  to  the  feed  point.  We  further  note  [I^]  = [I^] 

3 2 

and  [I  ] = [I  ] because  the  wire  geometry  is  identical  about 

its  vertical  axis.  Hence,  one  has  to  solve  for  only  unknown 

1 2 

current  vectors  [I  ] and  [I  ].  Also,  one  may  delete  the  third 
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and  fourths  rows  of  (4.3)  and  fold  the  resulting  matrix 
equation  of  two  by  four  into  a compact  partitioned  matrix 
equation  of  two  by  two: 


{(sj+c®]  + [cj+c^]} 

{[C^+C®]  + [C2+C2]} 


{[c^+cj]  + [C®+C®]} 
{[S2+C2]  + [C2+C2]} 


[I^] 

[Vi] 

[V2] 

(4.4) 


The  above  procedure  can  be  conveniently  applied  to  the  cases 
of  biconical  wire  antennas  placed  perpendicular  to  the  ground 
plane.  Fig.  4.3,  and  parallel  to  the  ground  plane.  Fig.  4.4. 

The  direction  and  symmetries  in  the  distribution  of  the  current 
in  the  corresponding  image  structures  are  shown  in  the  re- 
spective figures.  The  partitioned  matrix  equation  (4.1)  applies 
to  these  two  cases  of  biconical  antennas  perpendicular  and 
parallel  to  a ground  plane  with  a new  definition  for  the  matrix 
terms  , 


[cj]  = [C^]  + [C^]^  . (4.5) 


submatrices  obtained  for  a particular  element  and  its  corre- 
sponding image  element.  With  the  new  definition  for  the 
matrix  terms  in  (4.1),  step  by  step  matrices  can  be  folded 
and  (4.3)  can  be  solved  for  the  case  of  a biconical  antenna 
structure  placed  perpendicular  to  the  ground  plane.  Similarly, 
(4.4)  can  be  solved  for  the  case  of  biconical  antenna  structure 
placed  parallel  to  the  ground  plane. 
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TO  Ground  Plane 
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4.2. 


SAMPLE  RESULTS 


In  the  following  representative  sample  results,  viz.,  the 
distribution  of  the  current,  input  impedance,  and  far  field 
radiation  pattern,  are  presented  for  the  blconical  antenna 
structure  shown  in  Figs.  4. 1-4. 4. 

Fig.  4.1  shows  the  orientation  of  the  biconical  wire  antenna 

placed  vertically  over  a conducting  ground  plane.  The  structure 

is  assumed  to  be  in  the  xz-plane  with  the  feed  point  at  the 

origin  of  coordinates  on  the  conducting  ground  plane.  The 

geometry  of  the  wire  structure  is  specified  by  lengths,  h^  and 

h^,  and  radii,  a^  and  bottom  and  top  conical  angles, 

and  a2 » with  symmetry  maintained  for  the  other  two  wires.  The 

structure  becomes  a Bow-tie  antenna  [1]  with  02=180°.  The 

figure  also  shows  the  image  and  the  various  symmetries  of  the 

1 2 

antenna  currents.  Only  [I  ] and  [I  ] are  solved  for  this  geo- 
metry and,  in  Figs.  4.5  and  4.6,  are  shown  the  real  and  imaginary 
parts  of  the  current  distribution  on  a bow-tie  antenna  above  a 
ground  plane  at  five  different  frequencies.  The  current  is 
continuous  at  the  bend  and  goes  to  zero  at  the  upper  end.  The 
real  and  imaginary  parts  of  the  current  on  the  vertical  wire 
biconical  antenna  above  ground  is  shown  in  Fig.  4.7  at  three 
different  frequencies.  In  Fig.  4.8  the  input  impedance  on  a 
bow-tie  antenna  above  ground  is  plotted  on  a Smith  chart*for 
frequencies  from  400  MHz  to  2000  MHz.  The  input  impedance  makes 
one  complete  loop  in  this  frequency  range,  and  one  observes 
that  calculated  impedance  values  compare  quite  favorably  with 

*Smith  chart  data  are  normalized  with  respect  to  50  ohms. 
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Fig.  4.6.  Current  on  Bow-Tie  Antenna  above  Ground  ( a i=60°) 


9 3 


CURRENT  ( MA  ) 


• CALCULATED 
▲ MEASURED 


iG.  4.8.  Input  Impedance  of  Bow-Tie  Antenna  above 
Ground  ( Ct  ]^=60°) 


2a  = O.OI9 


thosi'  which  are  measured.  Further,  the  convergence  of  the 
input  impedance  of  the  bow-tie  antenna  as  a function  of  number 
of  coefficients  in  the  current  expansion  is  given  in  Fig.  4.9. 
Similarly,  in  Figs.  4.10,  4.11,  and  4.12,  are  shown  values  of 

input  impedance  of  other  bow-tie  antennas  above  ground  com- 
pared with  those  of  the  measured  results.  The  input  impedance 
of  a vertical  biconical  antenna  above  ground  is  indicated  in 
Fig.  4.13.  We  note  that  the  distribution  of  the  current  on  the 
vertical  wire  biconical  antenna  above  ground  and  the  wire  bi- 
conical antenna  in  free  space  is  identical,  but,  for  the  wire 
biconical  antenna  in  free  space,  the  input  impedance  is  twice 
that  shown  for  the  case  of  a vertical  biconical  antenna  above 
a ground  plane. 

For  convenience  and  ease  of  reading,  data  on  Smith  charts 
of  Figs.  4.8,  4.10-4.13,  4.15  and  4.17  are  tabulated  in 
Tables  III-IX. 

The  wire  biconical  antenna  in  free  space  is  now  placed 

above  a ground  plane.  Fig.  4.3,  with  the  plane  of  the  wires 

perpendicular  to  the  ground  plane.  The  feed  point  is  at  a 

distance  d above  the  ground  plane,  and  the  figure  also  shows 

the  image  structure  and  the  current  symmetries.  The  unknown 
12  3 4 

currents  [1  ],  [I  ],  (I  ] and  [I  ] are  determined  for  this 

geometry  and,  in  Figs.  4.14  and  4.15,  data  for  this  antenna 
are  given.  Similarly,  when  the  structure  is  place  above  and 
parallel  to  the  ground  plane  (Fig.  4.4)  at  a distance  d from 
the  ground,  the  corresponding  distribution  of  the  current  and 
the  input  impedance  are  as  shown  in  Figs.  4.16-4.17. 


A ;; .1  m |)  1 o far  field  pattern  of 


a wire  b i c o n i c a 1 ante  n n a 


in  free  space  is  given  in  Kig.  4.18. 


table  III 


Input  Impedance  of  Bow-Tie  Antenna 
above  Ground 
(Data  of  Fig.  4.8) 


Frequency 

Calculated 

Impedance 

Measured 

Impedance 

(.MHz  ) 

(Ohms ) 

(Ohms ) 

400 

9.8 

- 

j 

136.1 

12.4 

j 

00 

500 

12.3 

- 

i 

75.6 

12.6 

- 

j 

80. 9 

600 

14.  6 

- 

j 

31 . 2 

16.2 

- 

j 

33.2 

700 

24 . 0 

+ 

j 

9.0 

23.2 

+ 

j 

6 . 4 

750 

30. 5 

+ 

3 

29.5 

800 

39. 5 

+ 

j 

47.8 

38. 5 

+ 

i 

46  . 5 

900 

65.0 

+ 

,1 

95.0 

61 . 8 

+ 

j 

97.5 

1000 

108.5 

+ 

.1 

140.5 

97.5 

+ 

j 

140.0 

1 100 

180.1 

+ 

j 

192.0 

160.0 

+ 

j 

200.0 

1200 

35  5.2 

+ 

j 

186.3 

315.0 

+ 

j 

180. 5 

L 300 

525.7 

+ 

,1 

5 . 2 

510.0 

+ 

.1 

10.5 

1400 

390. 3 

- 

.1 

200. 1 

401.5 

_ 

,1 

195.0 

1500 

261.7 

- 

i 

275.  3 

280.0 

- 

j 

2 5 5.0 

1600 

105.6 

- 

1 

210.9 

155.5 

- 

j 

205.0 

I 700 

47.5 

- 

j 

155.2 

51.5 

- 

j 

165.0 

1800 

27.6 

- 

j 

110.4 

29  . 5 

- 

j 

115.0 

1900 

15.5 

- 

j 

61.5 

14.0 

- 

j 

65.5 

2000 

12.2 

- 

j 

15.5 

12.5 

- 

,j 

17.5 

2a  * 0.019 
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Convergence  of  the  Input  Impedance  of  Bow  Tie 
Antenna  above  Ground  ( Ct  i=60°) 
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Input  Impedance  of  Bow-Tie  Antenna  above 
Ground  ( O 2^=45°) 


BfcSiSTAHCt  COMPON^ 
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Input  Impedance  of  Bow-Tie  Antenna  above 
Ground  ( CL  i=30°) 


TABLE  IV 


Input  Impedance  of  Bow-Tie  Antenna 
above  Ground 
(Data  of  Fig.  4.10) 


Frequency 
(MHz  ) 

Calculated  Impedance 
(Ohms ) 

Measured  Impedance 
(Ohms ) 

600 

16.4 

- 

j 

39.4 

17.0 

_ 

j 

38.0 

700 

25.9 

+ 

j 

4.  2 

25.0 

+ 

j 

4 . 5 

750 

32  . 5 

+ 

j 

21.2 

33.5 

+ 

j 

24 . 0 

800 

40 . 5 

+ 

j 

42  . 2 

39.5 

+ 

j 

42.0 

900 

62.6 

+ 

i 

80. 3 

67 . 5 

+ 

j 

85.0 

1000 

105. 1 

+ 

j 

127.5 

113.5 

+ 

j 

120.0 

i 

1 

) TABLE  V 

Input  Impedance  of  Bow-Tie  Antenna 
above  Ground 
(Data  of  Fig.  4.11) 

1 


Frequency 

Calculated 

Impedance 

Measured 

Impedance 

i 

(MHz) 

(Ohms ) 

(Ohms ) 

1 

600 

15.  1 

- 

j 

61 . 3 

15.5 

_ 

j 

57.5 

» 

V 

700 

23.2 

- 

3 

18.6 

25.0 

- 

j 

14.5 

750 

28.6 

+ 

,1 

1.6 

32.0 

+ 

i 

6.2 

800 

35 . 5 

+ 

j 

21.  1 

39.0 

+ 

j 

27.5 

900 

53.  7 

+ 

j 

59  . 5 

47 . 5 

+ 

j 

52 . 5 

1 

1000 

81.4 

+ 

j 

98.6 

70.0 

+ 

j 

95 . 0 

s£SistA»*c£  ccvPCseiwT| 


CALCULATED 


MEASURED 


Fig.  ^1.12.  Input  Impedance  of  Bow-Tie  Antenna 
Ground  ( a i=15°) 


CALCULATED 


Fig.  4.13.  Input  Impedance  of  Vertical  Biconical 
Antenna  above  Ground  ( Ct  i=  q o=90°) 


TABLE  VI 


Input  Impedance  of  Bow-Tie  Antenna 
above  Ground 
(Data  of  Fig.  4.12) 


Frequency 
(MHz  ) 


Calculated  Impedance 
(Ohms ) 


Measured  Impedance 
(Ohms ) 


600 

700 

750 

800 

900 

1000 


15.7 

- 

j 

74  . 2 

24.2 

- 

i 

27.6 

29.5 

- 

j 

6.3 

36  . 3 

+ 

j 

15.0 

54.  2 

+ 

j 

51.5 

80. 5 

+ 

i 

99.4 

15.5 

- 

j 

82 . 5 

22.5 

- 

j 

35.0 

27  . 5 

- 

1 

12.0 

37.  5 

+ 

j 

7 . 5 

50 . 0 

+ 

j 

47.5 

85 . 5 

+ 

j 

87 . 5 

TABLE  VII 

Input  Impedance  of  Vertical  Biconical 
Antenna  above  Ground 
(Data  of  Fig.  4.13) 


Frequency 

(MHz) 

Calculated  Impedai 
(Ohms ) 

600 

33 . 9 

+ 

j 

40  . 7 

700 

68.6 

+ 

j 

98. 2 

750 

100.0 

+ 

j 

130.2 

800 

151.5 

+ 

j 

161.2 

900 

352  . 7 

+ 

j 

145.4 

1000 

411  . 5 

- 

j 

161 . 3 
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Fig.  ^.15.  Input  Impedance  of  Wire  Bi coni  cal  Antenna 

ABOVE  AND  PERPENDICULAR  TO  GROUND  PlANE 
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Fig.  4.17.  Input  Impedance  of  Wire  Bi coni  cal  Antenna 

ABOVE  AND  PARALLEL  TO  GROUND  PlANE 
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TABLE  VIII 


Input  Impedance  of  Wire  Biconical  Antenna 
above  and  Perpendicular  to  Ground  Plane 
(Data  of  Fig.  4.15) 


Frequency 
(MHz  ) 


Calculated  Impedance 
(Ohms ) 


600 

18.  9 

- 

j 

44  . 7 

700 

39.  7 

+ 

,1 

46.8 

750 

56.4 

+ 

1 

93.9 

800 

79.  9 

+ 

j 

142.9 

900 

157.9 

+ 

j 

244 . 4 

1000 

306.  3 

+ 

j 

335.9 

TABLE  IX 


Input  Impedance  of  Wire  Biconical  Antenna 
above  and  Parallel  to  Ground  Plane 
(Data  of  Fig.  4.17) 


F requency 
(MHz) 


Calculated  Impedance 
(Ohms ) 


600 

19.  5 

- 

j 

45.2 

700 

41.1 

+ 

j 

46 . 7 

750 

58  . 7 

+ 

3 

88.4 

800 

83.4 

+ 

j 

142.7 

900 

151.1 

+ 

j 

243.6 

1000 

319.5 

+ 

j 

324.6 

no 
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APPENDIX  1 


INTEGRALS  INVOLVING  PIECEWISE  LINEAR 
AND  PIECEWISE  SINUSOIDAL  FUNCTIONS 


In  solving  equations  by  numerical  techniques  incorporating 
piecewise  linear  and  piecewise  sinusoidal  testing,  one  encoun- 
ters the  need  to  evaluate  Integrals  involving  the  piecewise 
linear  and  piecewise  sinusoidal  "pulses."  In  this  appendix 
these  integrals  are  tabulated. 

I-l  PIECEWISE  LINEAR 


The  piecewise  linear  function  A (z)  is  defined  bv 

m ■' 


= \ (A-|z-z^|)  . ze(z_^_^,z^_^^) 


(1.1) 


w li  e r e 


A = (z  -z  ) . 
m m-1 


(1.2) 


is  depicted  in  Fig.  I-l. 


P iecewise  Linear  Testing  (Against  First  Derivative)  For  a 
function  f(z)  having  a first  derivative,  one  can  readily  show 
that  the  following  are  true. 
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(1.8) 


J = Ar"m+1>  " + ^("m-l> 


z = z 


m-  1 


1-2  PIECEWISE  SINUSOIDAL 


The  piecewise  sinusoidal  function  is  defined  by 


sin  k(A-| z-z  I ) 
Af<.) 


s in  kA 


and  is  illustrated  in  Fig.  1-2. 

Piecewise  Sinusoidal  Testing  (Against  First  Derivative; 
f(z)  possessing  a first  derivative,  one  can  derive  the 


ill 

/ 


4rf (z)A®(z)dz 
a z m 


f (z  ) 
m 


- k 


ill 

/ 


cos  k (z-z  , ) 
f(^)  dz 


s in  kA 


z-z 


m-  1 


z*z 


m-  1 


'm+  1 


I 


'm+l 


-f(z  ) + k 
m 


/ 


f (z) 


m+l 


sin  kA 


-d  z 


(1.9) 


For 

f ol 1 owi n g 


(I ■ 10) 


(1.11) 


' m+  1 


I 


^f(z)A®  (z)dz 
d z m 


sin  kA 


z = z 


m- 1 


/ 


f(z)cos  k(z-z  -)dz 
m—  1 


L 2 = Z 


m-1 


'm+1 


/ 


+ I f(z)cos  k(z  , -z)dz 
m+i. 


z = z 


(1.12) 


Piecewise  Sinusoidal  Testing  (Against  Second  Derivative) 

Finally,  for  f(z)  possessing  a second  derivative,  one  can  estab- 
lish the  relationships  which  follow. 


ni 

/ 


z = z 


„f(z)A®(z)dz 
dz^ 

m-  1 


- k‘ 


U1 

/ 


f (z)A“ (z)dz 
m 


(1. 13) 


'm-1 


/; 


m-f  1 
2 


— , f(z)A  (z)dz 

j ^ m 

dz 


^ c / ^ ,coskA-,  . k 

dz  sin  kA  sin  kA^^^m-Fl^ 


'm-Fl 


/■ 


z = z 


k^  J f(z)A®(z)dz 
m 


(1.14) 


m 
1 1 7 


Ei.V-.i-A.'. ’J* 


= HTTkA  ^ cos  kAf(z^)  + 


- k^  / f (z)A® (z)dz 

I m 


1-3  APPROXIMATIONS 


The  following  approximations  are  good  whenever  f(z)  is 
smooth  over  ^ ~ ) and  are  useful  in  applications  of 

piecewise  linear  and  piecewise  sinusoidal  testing. 


f (z)A^  (z)dz  = |-f  (z  ) 
m z m 


f (z)Af;’®  (z)dz  = |f(z  ) 
m z m 


f (z)A^ >®(z)dz  = Af(z  ) 

ID  m 


In  the  above  expressions  A 


represents  either  (I.l)  or  (1.9). 


APPENDIX  II 


ELECTRIC  SCALAR  POTENTIAL  DUE  TO  A RING  CHARGE 


In  this  appendix  is  presented  the  calculation  of  electric 
scalar  potential  due  to  a ring  of  uniform,  discrete  charge. 

An  expression  for  this  potential  is  often  needed  in  moment 
method  analyses  of  wire  scatterers  or  antennas  when  the  wire 
current  is  approximated  by  pulses.  A circumferentially  uniform 
surface  current  on  a cylinder  "deposits"  a discrete  ring  of 
charge  at  any  point  where  the  current  is  discontinuous.  If 
the  total  axial  current  on  the  cylinder,  assumed  to  be  coaxial 
with  the  z axis,  is  I^(z) , then  the  total  quantity  of  charge  Q 
in  the  ring  of  charge  is  given  by 


8. im 
alO  L 


(11.  1) 


where  z^  is  the  point  at  which  1^  is  discontinuous.  This  charge 
Q is  uniformly  distributed  over  the  ring  which  can  be  looked 
upon  as  a circular  line  charge  of  constant  linear  density. 


= _2_ 
2TTa 


(II. 2) 


where  a is  the  cylinder  radius. 

Knowing  the  density  of  the  ring  charge,  one  next  deter- 
mines the  scalar  potential  due  to  the  ring  charge.  If  the  ring 
charge  is  located  at  z=^  in  a cylindrical  coordinate  system. 
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with  the  plane  of  the  ring  parallel  to  the  xy  plane  as  depicted 
in  Fig.  I 1-1,  the  potential  <}>  at  a point  (r  ,<p  ,z)  is 


I 


I 


i 

1 


f 


I 


4>(r,$  ,z) 


TI 


4>  ' =-TT 


(II . 3) 


where  R is  defined  in  Fig.  II-  1.  Since  the  charge  is  constant 
in  <)>»so  is  0;  thus,  <I>  reduces  to 

$(r,z)  = (z-^;  , r , a)  (I  I.  A) 

where 


g (C  , r , a)  = 


^ / 


u 

2 2 2 ^ 
-jk[^  +r  +a  “2ra  cost})*] 

® d(D' 

[ 5^  + r^+a^-2ra  cos4>  ' ] 

= -7T 


(II. 5) 


A special  case  of  interest  requires  the  determination  of 
0(a,z)  in  which  case  g reduces  to  the  usual  kernel  (exact)  K(0 
used  in  thin-wire  calculations:  g(^,a,a)  = K(^).  Therefore, 

in  view  of  (II. 1),  the  potential  at  (a,z)  due  to  a current  dis- 
continuity at  z^  is 


Lastly,  it  is  worth  noting  that,  for  a<<A  and  |r-r'|>>a, 
the  potential  at  a point  r due  to  a ring  charge  at  r'  can  be 
written  simply  in  the  form. 


$ (a , z ) 


j^;^K(z-z^)Um 


rvin 


.J 
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-iklF-?' 


Him 

aiO 


. (11.7) 


again  with  the  ring  charge  resulting  from  a discontinuity  in 
^7.  2^. 

Singularity  of  Potential  Due  to  Ring  Charge  The  potential  $ 
due  to  the  ring  charge  is  singular  in  the  limit  as  the  point 
of  observation  approaches  the  ring.  Knowledge  of  the  behavior 
of  the  potential  on  a cylindrical  surface  due  to  a ring  charge 
at  the  end  of  the  cylinder  (Fig.  11-2)  is  useful  in  wire  anal- 
yses. Thus  we  investigate  the  manner  in  which  $(a,z)  becomes 
unbounded  as  Since  is  a constant  and  g(d,a,a)  = K(d), 

we  determine  the  limit  of  K(d)  as  d approaches  zero.  K(d)  can 
be  written 


K(d) 


7T 


TT 


2¥  /s'"*'- 


whe  re 


D = [ d^  + 2a^ (1-COS0  ' ) 


(II. 8) 


(II.  9) 


The  singularity  of  K(d)  is  due  to  the  first  integral  above 
and  the  second  is  well  behaved  so,  being  interested  in  the 
singularity,  we  concentrate  attention  on  the  first.  For  a 
positive  number  6<n, 
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and,  as  before, we  discard  bounded  quantities  and  retain  only 
the  integral  over  (-6,6),  which,  for  any  nonzero  6,  contains 
the  singularity.  Next  we  specify  that  6 be  fixed  but  suffi- 
ciently small  to  justify  replacing  2(l-cos<})')  in  (II. 9)  by 
2 

(J) ' . This  enables  us  to  reduce  the  integral  over  (-6,6) 
containing  the  singularity  to 

6 

4>'=-6  ^ 

which  can  be  readily  integrated  to  obtain 
In  the  limit  as  d-^0  , the  above  behaves  as 


TTa 


2 i m 
d'*’  0 


Hence  , 


K(d)  * 

d-^0 


1 

TTa 


c /2a. 
in(-^) 


(II.  10) 


and  the  limiting  behavior  of  the  potential  is 


1 


$ 


d >-0  ATI  ea 


(I  I . J J ) 


or,  as  a function  of  the  current  at  the  end  of  the  cylinder, 

2a 

0 \ £n(^)  . (11.12) 

d"^0  jAtt  cjca 


appendix  111 


W1  RF.-TO-WI  RE  DISTANCE  AND  LOCATION  SCHEME 


In  the  interest  of  flexibility  in  a wire  analysis,  one 
wishes  to  have  available  a general  scheme  for  specifying  the 
location  of  an  arbitrary  wire  as  well  as  the  geometric  pro- 
perties of  the  wire  wliich  are  pertinent  to  the  calculations 
to  be  performed.  Also,  one  must  be  able  to  specify  ir  a gen- 
eral manner  the  distance  from  a point  on  one  wire  to  that  on 
another.  A scheme  for  accomplishing  the  above  is  outlined  in 
this  appendix. 

As  suggested  in  Fig.  II  I-  1,  the  locations  of  the  lower  and 
upper  end  points  of  the  axis  of  a given  wire  are  denoted  in 

n ■ 1-  ./UUU, 

Cartesian  coordinates  by  (x  ,y  ,z  ) and  (x  ,y  ,z  ),  respective!; 
From  these  points,  one  determines  the  wire  length L,  the  direc- 
tion cosines  of  the  w'ire  cosa,  cosB,  cosy,  the  unit  vector 
along  the  wire  axis  u,  the  location  of  the  axis  center  r . 

Wire  Length: 


[(x^-x 


+ (y'^-yS"  + (z'^-z*'')^] 


(III.  1 ) 


Direction  Cosines 


cosa  = 


u «, 

X - X 


(II 1 . 2a) 
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(III.  2b) 


u 2 


u 2 

Z -7. 

cos  f = — 


Unit  Vector  (from  lower  to  upper  end): 


(I  I I . 2c) 


cosa  u + cosB 


+ cos  Y 


Location  of  Wire  Center: 


(III. 3) 


r 


c 


0 

X u 

X 


+ 


c ^ 

y % 


+ 


C /. 

z u 

z 


(III. 4) 


2 ^ L 

= X + — cosa  = 


1 / u 2, 

(x  +x  ) 


(III. 5a) 


c 2 L 1 u 2 

y = y + 2 cosB  = - (y^+y*^)  (III. 5b) 

c 2 L I u 2 

z - z + - COSY  = 2 (z  ) (III. 5c) 

Ihe  equation  of  a line  in  space  can  be  conveniently  formulated 
in  terms  of  the  location  of  a point  on  the  line  and  a parameter 
specifying  displacement  along  the  line  from  this  reference 
point  to  any  other  point.  With  the  wire  axis  center  as  refer- 
ence the  equation  of  the  wire  axis  is 


?(0  = + Cu 


(III. 6) 
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As  depicted  In  Fig.  III-l,  r(C)  locates  a point  on  the  wire 

— C 

axis  in  terms  of  the  fixed  reference  r and  the  displacement 
; u from  the  center  to  a general  point  on  this  axis. 

From  (III. 6),  one  may  readily  calculate  the  distance  from 
a general  point  r in  space. 


r=xu  +yu  +ZU 
X y z 


(III. 7) 


to  the  point  r'  on  the  wire  axis 


r - (r'^-i;  ' G) 


(x-x  ) - *cosa  u + 

^ X 


] “x  [(y-y^)  ■ ^'cosB  u^ 


■z  ) - i;'cosY 


Also,  one  may  determine  the  distance  D from  a point  on  the 

pq 

surface  of  the  p*"^  wire  to  a point  on  the  axis  of  the  wire 

(see  Fig.  3.2): 


D (? , C ' ) = 

pq  L 


2 i/”Ci,-^  \ / c 

a + (r+!;u)-(r- 

p p p q 


t h “ c 

In  (III. 9),  a is  the  radius  of  the  p wire,  r locates  its 
P P 


center,  and  u is  the  unit  vector  along  its  axis,  while  r 
P q 

locates  the  center  of  the  wire  whose  axis  is  directed 

along  u^.  The  parameters  and  denote  displacement  along 


p^^  and  q*”^  wires,  respectively,  from  their  respective  centers. 


APPENDIX  IV 


FAR  FIELD  COMPUTATION 


The  radiated  fields  of  a complex  coupled-wire  antenna 
structure  can  readily  be  calculated  with  far-field  conditions 
for  a known  current  distribution  using  either  reciprocity  or 
a vector  potential  approach.  In  this  analysis,  the  vector 
potential  approach  is  used  to  formulate  the  electric  field  in 
terms  of  current  distributions  computed  from  the  numerical  pro- 
cedure outlined  in  this  report.  As  would  be  expected,  this 
far-field  formulation  reduces  to  the  computation  of  an  equi- 
valent array  factor  and  pattern  factor  related  to  the  subdomain 
basis  set  used  to  represent  the  antenna  current  distribution. 

For  simplicity,  a single  wire,  as  the  p*"*^  wire  shown  in 
Fig.  IV- 1 , is  used  to  illustrate  the  formulation  which  can 
subsequently  be  extended  to  the  multiple-wire  case  required  to 
fabricate  the  wire  biconical  antenna.  From  standard  formulas, 
such  as  (3.3a),  for  total  vector  potential  in  the  direction  of 
the  current,  i.e.  the  wire. 


A (r) 
P 


0 


jk  r-r ' 

P 

d^  ' 

P 


(IV. 1) 


where 


1 30 


f-fpl  = (|r 


+ \r'l^  - 2r^r')^ 
P P 


and  Lp  is  the  total  length  of  the  wire,  as  defined  in 

Chapter  3 and  Fig.  IV. 1.  Now,  for  the  far-field  conditions. 


r large,  |r'|  bounded,  and  ■ 1 , (IV. 1)  reduces  to 

^ 1 r I 


A (?)  = 


P e 
47T 


f 

l?l  J 


j k I r ' I cos  $ 

ip(C')e  P ^dC'  (IV. 2) 


where  $ is  the  angle  between  r and  u . <P  can  be  expressed 

p p p 

in  terms  of  the  direction  cosines  related  to  r and  the  p*"^ 
wire  as 


cos<f  = cosa  cosa  + cosB  cosB  + cosa  cosa  . (IV. 3) 

P P ^ p r p r 

If  the  distance  |r'|  is  expressed  in  terms  of  and  the  dis- 

P p 

tance  along  the  wire  from  the  lower  and  to  the  upper  end,  C', 
the  vector  potential  becomes 


p 


h e 


-jk(|  ?|  -)  ?pj  cos'Jp)  ^ P 


/ 


j kC ' cos$ 

i^(? ')e  P d;' . 


(IV. 4) 


In  this  case,  r^  is  the  vector  from  the  origin  of  the  coordinate 

t h 2. 

system  to  the  lower  end  of  the  p wire  and  $ is  the  anele 

. P 

_ £ £ 

between  r and  r^  . <}>^  can  also  be  represented  in  terms  of  the 

related  direction  cosines  as 
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1 I I 

cosi  = cosa  cos  a + cosS  cos3  + cosu'cosot 
P P r P r P r 


(IV. 5) 


Kquation  (IV. 4)  is  simply  an  equation  for  the  total  vector 
potential  in  the  far-field  for  an  arbitrarily  placed  wire  with 
a general  current  distribution.  However,  if  the  current  dis- 
tribution is  expressed  in  terms  of  a linear  combination  of 
subdomain  elements  of  the  form  (2.6a)  and  as  shown  in  Fig.  2.3, 
then 


p ni  m 


(IV. 6) 


where  N is  total  number  of  subdomain  elements.  With  the  sub- 
stitution of  (IV. 6)  into  (IV. 4)  and  a subsequent  translation 
of  axis,  the  vector  potential  simplifies  to 


-jk(|  r|  - 


A (r) 
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r^  I cos<I>^ ) 
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j kmA  cos<I> 
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A (u)e 


j u cos  if 


du 


(IV. 7) 


-A 


where  A is  L^/(N-1),  N is  the  total  number  of  segments  on  the 


wire,  u is  distance  measured  relative  to  the  center  of  an  ele- 


£ £ 

ment  of  the  subdomain  basis  set  {A  },  and  A is  an  subdomain 

m 


element  centered  about  zero  on  the  u axis.  Note  that  at  each 
end  of  the  wire  the  integral  is  valid  only  for  A,  not  2A,  and 
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i 


U -1 


In  the  far-field  of  the  wire,  the  electric  field  in  the 
’ ‘^i’^^ction  of  the  wire  can  be  approximated  as 


r 


However,  the  6 and  C’  components  of  the  electric  field  in  respect 
to  the  origin  of  the  coordinate  system  are  the  qualities  of 
interest,  therefore,  the  projected  components  of  the  E-field 


Eg(f)  . 


(IV. 10a) 


E,(r)  = -JwA  (r)u  *u 
4>  p " p (; 


(IV. 10b) 


where  the  unit  vectors,  Uq  and  , can  be  represented  in  terms 
of  the  direction  cosines  related  to  r for  computation. 

A subroutine  PATRN  has  been  developed  for  the  far-field 
computation  of  the  total  0 and  (p  components  of  the  electric 
field  from  Eqs.  (IV. 8)  and  (IV. 10)  using  the  known  current  dis- 
tribution determined  for  the  geometry  defined  in  the  numerical 
procedure  presented  in  this  report  for  a single  wire.  The 
subroutine  computes  these  fields  for  the  piecewise  linear  sub- 
domain  basis  set  where  the  full-element  pattern  factor  is  of 
the  well-known  form 


p = A sin  (x/2) 
f (x/2) 


(IV. 11a) 


the  half-element  pattern  factor  is  of  the  form 


1 1 e+J* 

A - + J-  - ® 

a 2 - 2 

XXX 


(IV. lib) 


whe  re 


__=_======^  - ^ 

X = kAcosO 

P 

where  the  negative  sign  corresponds  to  the  upper  elements  and 
the  positive  sign  corresponds  to  the  lower  element. 

These  results  are  extended  to  the  case  of  multiple  wires 
encountered  in  the  wire  biconical  antenna  by  linear  super- 
position of  the  total  electric  field  as  indicated  in  (3.1). 


136 


